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NAVAL ARCHITECTURE 


The current and future possibilities 
to adapt the R-JH3 method 
to ship power prediction practice 


Henryk Jarzyna 
Institute of Fluid-Flow Machinery, 
Polish Academy of Sciences, Gdansk 


ABSTRACT 


In the paper the current and future possibilities of the R-JH3 method are discussed. The R- 
-JH3 method (presented in detail in [27]) was elaborated for the ship power prediction on 
the basis of a new definition of the open-water propeller fully equivalent to the behind-the- 
-hull propeller. The R-JH3 method can be used to determine the mean effective velocity and 
the mean effective pitch coefficient of the behind propeller. The R-JH3 method may be very 
useful for scale-effect investigations and scale-effect designation when results of model 
self-propulsion test have to be transformed to the ship scale. 


Keywords : ship hydromechanics, theory of propellers, ship hull-propeller interaction 


PECULIARITY OF THE ITTC-78 
POWER PREDICTION METHOD 
AND NECESSITY OF CHANGE 
OF THE THIRD GROUP OF PROCEDURES 


The ITTC-78 ship power prediction method 
can be divided into three groups of procedures : 


> the 1* group of procedures dealing with the ship resistance 
Rts 

> the 2™ group of procedures concerning the propeller thrust 
Tps necessary to ensure the ship velocity Vs 

> the 3™ group of procedures for prediction of the propeller 
torque Qpgs and the propeller revolution number ng. 


Only the a group of procedures 
is proposed to be changed. 


In the 1 group of procedures the ship model resistance 
is determined and its experimentally obtained value is trans- 
formed to the full scale. The scale effect is the main problem of 
the transformation. 


In the 2"¢ group of procedures the change of the ship 
model resistance due to the action of the behind propeller is 
determined in the form of the suction coefficient. The value of 
the model suction coefficient is further supposed to be the same 
as that in the full scale. 


In the 3"4 group of procedures being the subject of chan- 
ges, the ITTC-78 method introduces the concept of the equiva- 
lent open screw propeller. The concept makes it possible to 
determine — in model scale — the effective mean velocity and 
the advance mean coefficient. The possibility can be realized 
when the model thrust coefficients of the behind propeller and 
the open propeller are equal to each other. Their identity gives 
the value of the mean advance coefficient Jry and the effecti- 
ve mean velocity VtM = JTM ` om: Dm. 


The torque coefficients of the two propellers for J = Jt) 
are different. The ratio of the two torque coefficients values : 


Keom (J ™ ) 
Kosm 


RM 


is known as the relative rotative efficiency. 


The effective mean velocity Vrm and the open propeller 
(SOM) characteristics Krom = f(J) and Kgom = f(J) are trans- 
formed from model scale to full scale according to the ITTC- 
-78 transformation formulae. This way the mean effective velo- 
city Vrs is received. The hydrodynamic characteristics Ktos = 
= f(J) and Kaos = f(J) are to be joined with a new open pro- 
peller in full scale. In model scale the identity of the thrust 
coefficients of the model behind propeller (SBM) and the equi- 
valent model open propeller (SOM) : Kppy = Ktom(J), is the 
criterion of equivalence of both model propellers. In full scale 
the identity of the thrust loading coefficients of the behind pro- 
peller (SBS) and the open propeller (SOS) is the equivalence 
criterion of these propellers : 


Tss __ Kros (J) 
2 n2 2 
Ps VrsDs J 
where : 

s z~“ thrust loading coefficient of the SBS 
PsVīsDs propeller 
Kros W) () - thrust loading coefficient of the SOS 

J propeller. 


From this equivalence criterion the advance coefficient 
J = Jrs can be determined and then the number of revolutions 
of the propeller, ng : 
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Vrs 


n g= 
JrsDs 
In the ITTC-78 method it is further necessary to assume 
that the value of the relative rotation efficiency in full scale, 
NRs, is equal to ngm. From the relationship : 


| Koom (Im) Koos (rs) 


RM 7 
Koem Kops 
the torque coefficient Kogs can be determined : 
K= Koos (J TS ) 
QBS 7 
NRM 


and hence the torque Qgs and the power Pps: 
2 5 
Qrgs = Kags “Ps “Ns “Dg 


Pps = 27 : Qgs ' ns 
can be calculated. 


One can note that the 34 group of procedures being in use in 
ship power determination has been designed on the basis of the 
equivalent open propeller definitions different in model and full 
scale. The geometry of the SOM propeller and the SBM propel- 
ler are identical. The a priori defined geometry leads to the 
hydrodynamic characteristics of such open propeller, Kpo = 
= f(J) and Kgom = f(J), and to the equivalence criterion limi- 
ted to the thrust coefficients only : 


Ktpm = Kromi) 

From the criterion the advance coefficient J = Jyry results. 
Hence the effective mean velocity Vrm = Jrm ` nm: Dm can 
be determined. In full scale the hydrodynamic characteristics 
Kros = f(J) and Koos = f(J) derived form the model relation- 
ships : Krom and Kgom with the help of the ITTC-78 trans- 
formation formulae, define the open propeller (SOS). 

The equivalence criterion of the SOS propeller and the SBS 
propeller is different from that in model scale and is given in 
the form of equality of the thrust loading coefficients of the 
SBS and SOS propeller. 


Some controversial elements of the 3° group 
of procedures of the ITTC-78 method can be indicated. 


O The open screw propeller is given a priori being different 
in model and full scale : 


+ in model scale : (SOM) = (SBM) 
+ in full scale : the SOS propeller is described in the form 
of the characteristics Kaos(J) and Kyos(J). 


OQ The equivalence criterion of both propellers (the behind 
propeller and open one) is different in model scale (the equa- 
lity of thrust coefficients) and in full scale (the equality of 
thrust loading coefficients). 


Q The requirement that only the thrust coefficients of both 
propellers are to be equal to each other in model scale, re- 
sults in that the torque coefficients of the propellers are dif- 
ferent. The ratio of the coefficients, i.e. the relative rotative 
efficiency, is the same — according to the assumption — in 
model and full scale. Due to such assumption only it is 
possible to determine the torque Qgs and the power Pps by 
means of in the ITTC-78 method. 


O The condition of equality of the thrust loading coefficients 
of both full-scale propellers is based on the following as- 
sumptions : 
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+ the transformation of Vrm into Vrs is known 
+ the transformations of Krom and Koom into Kros 
and Koos are known. 


The transformations can be verified 
by using indirect methods only. 


The controversial elements in question 
can be eliminated, in principle, in two ways : 


# by introducing the generalization of the equivalent open 
propeller definition 

# or by eliminating the equivalent open propeller definition. 
In this case the numerical determination of the ship power 
can be taken into account in the full scale directly. 


The first way has been applied in the Jarzyna’s papers [20, 
21, 22, 27]. The second one is wholly connected with investi- 
gations into the computer hydrodynamics problems (numeri- 
cal basin). It is a very promising direction of scientific activity 
though its practical application is not yet possible. 

Even though this research direction will bring full success, 
the results of investigations into the equivalent open water pro- 
peller will remain of a great value in determining the effective 
mean values connected with the behind propeller problems. 

As an example of such values can serve the effective mean 
velocity or the effective mean pitch coefficient of the behind 
propeller. The only way to determine the mean effective valu- 
es connected with the behind propeller is to make use of the 
definition of the equivalent open water propeller. 


So the generalization of the equivalent open propeller 
definition can satisfy the following different conditions : 


# to make the realization of new concepts of the 3" group of 
procedures of the ITTC-78 method possible 

* to make it possible to form a compact criterion of determi- 
nation of mean effective values which are connected with 
screw propeller’s global characteristics 

# to serve as a valuable tool in scale effect investigations espe- 
cially when the propeller thrust full-scale measurements 
are possible in practice. 


THE ESSENCE OF THE R-JH3 METHOD 


Definition of the open water propeller 
equivalent to the behind propeller 


For given values (Krg , Kap) of the behind propeller and 
for the one — parameter family of open — water screw propel- 
lers with the pitch coefficient (P/D) being the family parameter 
(of the same family geometry as that of the behind propeller, 
except of the propeller pitch) the equation set : 


P 
() = Kop = Kyo {5 


P 
Q) > Ky=K 5] 
QB QOj| `> 
D jJ; 


has one and only one solution [J = Jro ; P/D = (P/D)e]. The 
parameter (P/D), identifies — out of the propeller family — one 
propeller SO[(P/D).] which is the wanted open propeller (SO) 
equivalent to the behind propeller (SB). The characteristics of 
the SO[(P/D).] are : 


-Pi 


rool, 


Identity of other global characteristics 
of both propellers 


The global hydrodynamic parameters of both the propel- 
lers: SB and SO[(P/D),] satisfy the following equations obtai- 
ned from the identity of thrust loading coefficients and the iden- 
tity of torque loading coefficients : 


5 - 
6) => = 
PVigD J 

P 
Q Ko) 
(4) => — 4. = +} + 
PVrgD J 


In order to obtain Vro from the equation (3) the value J = Jro 
resulting from the set of equations (1) , (2) is to be used. 


Thereafter the value of the number 
of propeller revolutions, n : 


can be received. 


(56) > n= 


The wanted value of the torque Qg results from the equa- 
tion (4) for J = Jro hence the delivered power Pp can be calcu- 
lated as : Pp = 27 - Qu: n. 


Next the thrust and torque coefficients 
Krs and Kog can be verified. 


THE PROBLEMS CONNECTED 
WITH THE USE OF THE R-JH3 
METHOD IN FULL SCALE 


The set of equations (1), (2), (3), (4) and (5), in presence of 
the feedback between the parameters Krgs , Kaps and the 
roots of the equation set, can be solved by means of the itera- 
tion method only. One can suggest two procedures to determi- 
ne the first approximation values of the Krgs and Kags. 


* The first procedure is based on the preliminary design of 
the behind propeller. To this end the following input para- 
meters should be given : 


Tgs - the propeller thrust sufficient 
to maintain the velocity Vs 

Vs - ship velocity 

Dg - screw propeller diameter 

V_/Vs(t/R) - the design velocity field. 


The thrust Tps results from the first and second procedure 
group of the ITTC-78 method. The parameters Vs and Ds 
are given values. The design velocity field can be produced 
by using the Koronowicz’s procedure (computer program). 


* The second procedure of determination of Krgs and Kags 
is the O-JH3 method of ship power prediction. 


The O-JH3 method is below presented in detail because it 
can serve as an independent prognostic method parallel to its 
function in determining the first approximation values of KTgs 


and Kags for the R-JH3 method. The O-JH3 method based on 
the new definition of the equivalent open water propeller is 
very strongly connection with the ITTC-78 method and it can 
serve as a very useful element of correlation of both power 
prediction methods. 

It is worth mentioning that in the case when the paralell 
measurement results of thrust and torque, number of propeller 
revolutions and ship velocity are at one’s disposal, the R-JH3 
method will be a useful tool for verification of and for scale- 
-effect investigations on selected parameters applied in ship 
hydrodynamics. 


GENERAL REMARKS CONCERNING 
THE REALIZATION OF THE PROGNOSTIC 
METHODS : R-JH3 AND O-JH3 


> The realization of the methods R-JH3 and O-JH3 does not 
entail any additional experimental activity in comparison 
with that related to the ITTC-78 method. 


> In the 3 group of procedures the results of experimental 
model tests are processed by means of computer programs 
quite different from those used in the ITTC-78 method. The 
proposed methods do not use the open-propeller hydrody- 
namic characteristics experimentally obtained as in the case 
of the ITTC-78 method. In the new method the hydrodyna- 
mic characteristics of the one — parameter open — propeller 
family, and thereafter of the equivalent open propeller, are 
determined by using a computer program. The propeller 
family’s geometry is the same, except of the pitch coeffi- 
cient equal to that of the behind propeller. 


> The realization of the R-JH3 and O-JH3 methods, in the 
phase of their introduction to the prognostic practice, is 
carried out parallel to the calculation procedures of the 
ITTC-78 method. 


> Inthe ITTC-78 method the value of the parameter : 
Vrs 
Vs 

Vim 


Vm 


is forejudged because the mean effective velocity Vrs is 
determined from the model value Vry by using the ITTC- 
-78 transformation. 


Go = 


In the O-JH3 method the variability interval of the parame- 
ter G is assumed to be Gmin < G < Gmax, in which the value 
Goo is contained. The results of the O-JH3 method are given in 
function of the parameter G (G = Gog included). 


In the R-JH3 method the effective mean velocity Vrs re- 
sults directly from the equation (3) (due to the identity of thrust 
loading coefficients of the behind propeller and the equivalent 
open propeller) for J = Jros and P/D = (P/D)es being the roots 
of the set of equations (1) and (2). 


In the R-JH3 method the mean effective velocity Vros is 
directly dependent upon the input values (KTgs,Kosgs) to this 
equation set. 


THE REALIZATION OF THE R-JH3 
METHOD FOR GIVEN (Kg , Kog) VALUES 
OF THE BEHIND PROPELLER 


One can distinguish three cases 
when the (KrB,Kop) values can be determined univocally : 
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> Full-scale measuring results : Tgg, QBs, Vs, Ds - 


In the case when the ship measuring results (Tps , Qgs , ng, 
Vs) are available, then — for given Krgs, Kops — the equiva- 
lent open propeller, SOS[(P/D),,], can be received from the 
equations (1) , (2) with the roots [Jtas,(P/D),s]. And, Vros 
results from equation (3). Then the number of revolutions 


__Vros 
JtasDs 
and the torque Qpgs can be calculated from (4). 


Ng 


Both the values ng and Qgs should be comparable with tho- 
se measured on ship. This way a valuable correctness crite- 
rion of the measured values was created. The iteration me- 
thod can be used to bring the values to identity. 


On the other side the full-scale values : 


Vros , Kros = £ ([J,(P/D)es] 
Koos = f[J,(P/D)es] , (P/D)es 


can be compared wih those in model scale : 
Vtoem > Krom = f ([J,(P/D)em] 
Koom = f [J,(P/D)em] , (P/D)em 


The method R-JH3, used in parallel for both model and 
full scale, can be a very efficient tool in studying scale-effect. 


Ə Model-scale measuring results : Tgm > QBm > Vm > Ou - 


The R-JH3 method can be used in model scale because the 
values measured in self propulsion tests make it possible to 
determine Krgm and Kogm. The set of equations (1), (2), 
(3), (4) and (5) gives, in model scale, the following results : 


Krom L.@/D)em] 
Tro: P/D)em —> SOM[(P/D) 1] > Kooml,P/D)eml 


VTM >M ; QBM 
Two of the values, ny and Qpy , if related to the measured 
values, can be a measurement correctness criterion for the- 
se values. The iteration method can be used to bring the 
values to identity. 


> Calculated values : Tgs , Qgs , Vs , ng 
for a given final propeller (SBS). 


The R-JH3 method makes it possible to find — for a given 
behind propeller and given parameters (Kgs , Kags) — the 
following mean effective values : 


= the effective mean velocity Vro 
= the effective mean pitch coefficient (P/D),. 


THE R-JH3 METHOD USED 
TO POWER PREDICTION DIRECTLY 
IN FULL SCALE 


The difficulties in use of the R-JH3 method directly in full 
scale have been already mentioned. The difficulties are due to 
the feedback of the behind propeller and open propeller para- 
meters. Two possible ways of overcoming the difficulties were 
indicated. 


1* — which concerns the preliminary behind propeller design 
used as a source of the first approximation values of the 
thrust and torque coefficients : Kres , Kags. 

One can assume that the designed propeller is the optimum 
one from the efficiency point of view. For a given propeller 
diameter an optimum revolution number can be determi- 
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ned. The unknown design velocity field can be described 
by using different approximate methods, e.g. the results 
presented by T. Koronowicz [28]. In the referred papers the 
effective velocity field is numerically generated as the final 
result of interaction between the ship hull and the screw 
propeller. Although the investigations are still in progress 
their partial results can be used when only the first approxi- 
mation values are needed to start using the R-JH3 method. 


24 _ concerning the zero method O-JH3 used to produce the 
first approximation values of Krps , Kags. 


Stage 1 : the R-JH3 method is used in model scale. The 
parameters Jrom , (P/D)em and the equivalent open propeller, 
SOM[(P/D) em] , of the hydrodynamic characteristics Ktow 
and Kgom , are received. The effective mean velocity : 

VTOM = JTQM - np ` Dy can be determined. 


Stage 2 : the parameters VToM , Krom and Kagom are 
transformed to full scale by using the ITTC transformation 
formulae. 


Stage 3 : from the equations (3) , (4) and (5) of the R-JH3 
method applied to full scale the parameters Jros , ns and Qgs 
can be calculated. 


It should be mentioned that Vros can be determined — ac- 
cording to Griegson and Holtrope [5 , 7] — as a variable para- 
meter given by the relation : 


Vtos = G: VTOM 

where G is defined as : 

Vros 
G= ; 


Gmin < G < Gmax 


Vm 
In the case when the variable parameter G is introduced 
all results obtained in stage 3 are functions of this parameter. 


Thereafter the thrust and torque coefficients Kgs and Kogs 
are calculated in order to be used in the R-JH3 method as the 
first approximation of the input data. 

The values Qgs = f(G) and ng = f(G) 
can be used to predict the delivered power Pps: 
Pps = 27 : Qgs ` ns 
NOMENCLATURE 
G 


the ratio of nondimensional effective mean velocity of 
ship and model 


hr - relative rotative “efficiency” 

J - advance coefficient 

Ko - torque coefficient 

Kr - thrust coefficient 

nm» ng  - number of revolutions of model and ship propeller, 
respectively 

P - power (in general) 

Q - propeller torque 

SB - behind - the - hull (behind) screw propeller 

SB roa - modified behind screw propeller with constant 
effective pitch coefficient 

SBM - behind propeller in model scale 

SBS - behind propeller in full scale 

SO - open-water propeller 

SO[P/D)e]- open-water propeller with the given effectiv pitch 
coefficient 

SOM - open-water propeller in model scale 

SOS - open-water propeller in full scale 

T - propeller thrust 


Vs - ship velocity 
Vro - mean effective velocity from thrust and torque 
identity 


Indices 


B - behind quantity 

M - model quantity 

0 - open water quantity 
S - ship quantity 
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An algorithm for preliminary estimating 
hull structure mass and mass centre height 
of inland navigation ships 


Jan P. Michalski 
Gdańsk University of Technology 


ABSTRACT 


oe 


The paper presents an algorithm for preliminary calculation of mass and mass centre 
height of hull structure of inland navigation ships. It was elaborated basing on require- 
' | ments of the Rules for the Classification and Construction of Inland Navigation Ships of 
á Polish Register of Shipping, with application of a simplified method of estimating mass of 
hull plating stiffeners. The algorithm deals with the dimensioning of scantlings of structural 
members of classical ships intended for shipping dry cargo, and dry and liquid bulk cargo- 
es; however it does not cover ships of entirely different structural arrangement such as roll- 


-on-roll-off ships fitted with heavy decks, as the dimensioning of their scantlings is based on different models. 


Key words : preliminary ship design, ship structure, preliminary ship mass estimation 


INTRODUCTION 


In the ship design theory an important role is played by the 
methods intended for preliminary, approximate estimation of 
mass and mass centre coordinates of a designed ship. 

Knowledge of the parameters is necessary already in early 
designing phases to balance ship floatability and stability be- 
fore exact structural strength calculations are performed. In the 
subject-matter literature many similar methods dealing with 
sea-going ships are known e.g. [1, 2, 3, 4]. However only few 
methods concerning inland navigation ships have been publi- 
shed so far, e.g. [5, 6, 7]. 

Such methods are elaborated on the basis of statistical in- 
vestigations of mass of built ships or the rules of classification 
institutions, which provide analytical and tabulated relation- 
ships or diagrams. Applicability range of a particular method 
is usually limited to ships of a given functional type and hull 
structural arrangement; moreover it concerns a given structu- 
ral material and ship size range. 

For elaboration of the algorithm in question the rules for 
dimensioning hull structural elements, contained in the Rules 
of Polish Register of Shipping [9], were applied, as well as 
a method of accounting for mass of hull plating stiffeners [5] 
was used. The simplyfing assumptions applied in the conside- 
red mathematical model have been deemed allowable for ap- 
proximate estimating the total mass of ship hull structure. 

The elaborated calculation algorithm has been made easy 
for computer programming as it is intended for preparation of 
a computer software to perform series of simulation calcula- 
tions of ship hull structure mass to provide a basis for elabora- 
tion of a parametric method to predict mass characteristics of 
inland navigation ships. 
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An origin of the presented research are problems met du- 
ring the design work on realization of the Eureka - INCOWA- 
TRANS (Inland and Costal Water Transport System) project 
E!3065. One of the goals of the project is to develop a modern 
fleet for ecological inland waterways navigation. 


Formulation of the problem and its assumptions 
The algorithm concerns inland navigation ships intended 


for shipping both solid and liquid cargoes, whose construction 
corresponds with the midship section shown in Fig.1. 


JL 
| JE ES 


Fig. 1. Schematic diagrams of standard midship sections 
of inland navigation ships 


The presented standard midship sections are special cases 
of the generalized midship section shown in Fig.2. By respec- 
tive establishing the following dimensions : double bottom 
depth, double side breadth, hatch opening breadth and hatch 
coaming height as well as number of longitudinal bulkheads — 
— an appropriate standard midship section can be generated 
from the generalized one defined in the algorithm. 


Fig. 2. Schematic diagram of the generalized midship section 
of inland navigation ships 


The considered research problem consists in formulation 
of a mathematical model for analytical transforming the ship 
characteristics known at preliminary design stage, namely : 


% ship numerical parameters X = (L, B, H, T, Cp =) describing : 
ship main dimensions, hull form, structural material pro- 
perties 

% qualitative attributes P = (K, R) identifying : ship functio- 
nal type, its class and midship section topology 

¢ into values of the hull structure mass m and mass centre 
height z, , to be determined : 


m(x)=m/(L,B,T,H,Cg.p;,.--,P,) (1) 
z @)=2(L,B, T, Ce Pis Pa) 2) 


The hull structural elements should satisfy a number of 
strength requirements imposed by the classification institution 
rules. In the algorithm in question dimensioning the scantlings 
is carried out on the basis of overall and local strength criteria. 


ALGORITHM FOR DETERMINING 
HULL STRUCTURAL SCANTLINGS 


Qualitative identification attributes 
of hull structures 


Two functional types of ships, identified by means 
of the symbolic index K, are considered : 


> K= 1 dry cargo ships 
> K=0 liquid cargo ships. 


Three classes of ships determined by restricted areas of 
navigation and identified by the symbolic index R [10], are 
distinguished : 


* R= 1 -— navigation area on which waves of the height 
h, /;9 < 2.0 [m] can be encountered (e.g. Bay of Gdansk) 


* R= 2 -— navigation area on which waves of the height 
hy 49 <1.2 [m] can be encountered (e.g. Vistula Bay, Wło- 
cławek Bay) 


* R= 3 - navigation area on which waves of the height 
hy 49 £0.6 [m] can be encountered - rivers, channels and 
lakes, which are rated, by the administration, among the 
inland waterways of Poland. 


Numerical identifiers of hull structures 


Determination of ship hull scantlings and its structural desig- 
ning is associated with knowledge of ship hull form which — 


— at early design phases — is identified by : ship dimensions, 
hull space subdivision, hull form coefficients etc. In the algo- 
rithm in question the following numerical parameters identi- 
fying the hull were assumed : 

L - ship length [m] 

B - ship breadth [m] 

H - ship depth to upper deck [m] 

T - ship design draught [m] 

Cp - hull block coefficient [-]. 

Design loads 


The longitudinal hull bending moment Mpw is the sum of 
the still water bending moment Mr and wave bending moment 
Mw. In accordance with the requirements [9] : 


æ% the still water bending moments were determined as that 


not less than : 
M; =0.07B-H-L* [kNm] (3) 
æ for determination of the wave bending moment the follo- 
wing approximate formula was used : 


Mw = k,(R)-k, (T). AR, L) 


ehi (R):Cp B-L [kNm] 

where : h;(R) - design wave height. 
The correction coefficients kı(R) , k2(T) , Ay(R,L) in [9] 
are defined by means of tabulated discrete values, rather incon- 
venient in coding computer algorithms. Therefore the relation- 


ships were approximated by using the following analytical 
functions : 


h,(R)=h(R)=3-1.1R +0.1R?° (4) 
(R)= f£,(R)=0.81+0.095R —0.055R* (5) 


k 
k, (T)=f,(T)=2.2+0.014T? -1.85T°? (©) 
A,(R,L)= f;(R,L)=0.37-0.0115R + 
—0.0016L —0.00063R -L 
The design bending moment is the sum of the moments : 
Mc = Mew =M; +My [kNm] (8) 


Permissible stresses 


Ship hull structure exposed to external loads due to action 
of environment during building and operation of the ship, sho- 
uld provide an appropriate level of its resistance against : 


> loss of integrity of structural elements 
> excessive geometrical deformation of the structure 
> permanent change of dimensions of the structure. 


Scantlings of structural elements should be so selected as 
to fulfil the above mentioned criteria at given material proper- 
ties, that can be realized by an appropriately low level of inter- 
nal loads within the structure. 

The level of permissible stresses in the structure was defi- 
ned by means of the coefficient k < 1 related to the structural 
material yield stress op [MPa] . Hence the required strength 
modulus of midship section can be described by the following 
formula : 

Mc 


W, =10 
k 


[cmm] (9) 
. Op 

As the approximating formulas usually are of a nonstructu- 
ral character hence they are provided with appropriate units of 


measure. 


POLISH MARITIME RESEARCH, No 1/2005 9 


NAVAL ARCHITECTURE 


NAVAL ARCHITECTURE 


Description of midship section geometry 
is identified by the set of parameters : 


e - distance from the neutral axis [m] 

B, - breadth of cargo hatches [m] 

Ba» - breadth of inner bottom [m] 

H, - height of hatch coamings [m] 

Hg, - depth of double bottom [m] 

tẹ - thickness of outer bottom and side plating [cm] 
tą - deck plating thickness [cm] 

tg, - thickness of inner bottom [cm] 

tg; - thickness of inner side plating [cm] 

ty - thickness of longitudinal bulkhead plating [cm] 
S - midship cross-section area [m°] 

We - minimum midship section modulus [cm’m] 


Minimum thickness values of the elements taken into ac- 
count in the section modulus W, cannot be smaller than the 
minimum ones defined in [9], where a - frame spacing [m] : 


minimum thickness of deck plating : 
ty =a(0.01L+0.5) [cm] (10) 


minimum thickness of outer bottom and side plating : 


t, =a(0.01L+0.65) [cm] (11) 
minimum thickness of inner bottom plating : 
ta =a(0.01L+0.45) [cm] (12) 
thickness of inner side plating : 
ta =1.77a(0.082L+2.5) [em] (13) 
thickness of longitudinal bulkheads : 
tp =5avH +0.1 [cm] (14) 


The midship section modulus at minimum thickness 
of deck plating is determined as follows : 


2.4H? 


f()(142.12 


The midship section modulus at minimum thickness of 
bottom plating is determined by means of the following for- 


W, =t [cmm] (15) 


mula : 2 
min 2.4H 
a T 
f0(1+212] 
: 16 
£) B oe (16) 
1-f;() B-B, 


The value of the dimensioning section modulus Wc 
is the largest out of the following ones : 


min(W, )= max(W,, W, W3) (17) 


Mass of hull structure stiffeners 


The functions f4 and f5 which represent magnitude and di- 
stribution of mass of hull plating stiffeners, given in [5] by 
means of diagrams ( rather inconvenient for coding computer 
programs) are described by the relationships : 


r0=1 7? 2) 


H H 


(18) 


jie i independent variable 
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— discretely changeable parmeter. 


B ty, H 
f; () =f; E ca > He 
d 


— independent variable 


T| w 


(19) 


B t 
where: — $ „w, 


h ta 


oh — discretely changeable parameter. 


The relationships were approximated 
with the use of the analytical functions : 


0.95 P 0.82 
B-B, B 
H H 


ds| 0.05 + 


-107 
(20) 


fy () = 


1.68 
1-ds)| 0.145 + —— 
+(1-ds) +55, 

H 

where : ds = {0,1} 


If ds = 1 then inner side structures are taken 
into account in the hull section modulus. 


0.125 
Bao tan 
B-B,. tj (21) 
where : db = {0,1} 


If db = 1 then inner bottom structure is taken into account 
in the hull section modulus, and in the case of single bottom 
hull : f5(-) = 0.5. 


Hay 


f,(-)=db| 0.254 +0.55 ral 0.5(1—db) 


Mass and mass centre height 


The cross-section area corresponding with the midship sec- 
tion modulus Wc is determined by means of the formula con- 
sisting functions which account for stiffeners. In the case of 
a double bottom hull, the cross-section area of midship section 
is determined by means of the relationship : 


B-B W, H 
(22) 


H H 


B ta +t 
4+£,|—,—_—* [Bata [m 
(2 {i db *db 


In the case of a single bottom hull, the cross-section area 
of midship section is determined by the formula : 


Saai ae 2) (142.158), 
(23) 


H H 


{1+1342},o027H(,, +) [m°] 


In [5] the functions fe and f7 are defined by means 
of the diagrams expressing the relationship : 


soaa 2) 


H OH 


(24) 


—B, 


where : B - independent variable 


5 - discretely changeable parameter. 


B ta +t 
f ()= piate 


tap 


(25) 


where : z. independent variable 


ty. +t : 
oe discretely changeable parameter. 


db 


The relationships were approximated by using 
the following analytical functions : 


i138 + 0252 - us) 
fe()= B-B 10“ 
+(t-as{17 + 012) 


1.22- 0.125 È +3.14 to tte y 


tap 


=).375 B tastte 
H 


10° 
(27) 


f, () = 

tab 

The total mass of hull structure is expressed by the formula : 
m(X,B)= Pw S. P)L-fs(K, R)I 


where : 


(28) 


fg C) - a function for correcting mass 
distribution in hull end regions 


Pm [tm°] - structural material density. 


The function fg(-) defined in [5] by means of tabulated values, 
was approximated by using the analytical function : 


f,(K,R)=1.92+0.167K-0.225R (29) 
The height of hull mass centre Zg (+) 
is defined by the relationship : 
e=z,()=H-f5()=H- (30) 


«| db} 0.25+ ue 
Bab tab 


+0555 +0.5(1-db) 


where : the function f5(-) is defined by the formula (21). 


PRELIMINARY VERIFICATION 
OF THE ALGORITHM 


The preliminary verification of the algorithm was perfor- 
med by using — as an example — SINE 207 ship designed in the 
frame of the Eureka E!3065 INCOWATRANS project. The se- 
lected ship is intended for carrying general cargo and contai- 
ners : one tier in the hold and another on the deck hatch covers. 


Main parameters of the SINE 207 ship 


overal ship lenght Loa = 56.50 m 
ship lenght between perpendiculars Lgp = 55.30 m 


ship breadth B =9.00m 
ship depth to upper deck H =3.00m 
design draught Tp =1.00m 
deadweight at Tp Pp =280t 
scantling draught Ts =1.60m 
deadweight at Ts Ps =580t 
container capacity Cc = 18/36 TEU 
power N =620kW 
speed v =15km/h 


In the ship's design documentation, the outfitted ship mass 
equal to 144.5 t was estimated on the basis of the longitudinal 
distribution curve of ship mass. The abscissa of mass centre 
was estimated equal to 27.66 m, and mass centre height equal 
to 1.95 m. Because of lack of information on a navigation area 
assumed in the design of the ship, the verifying calculations 
were performed for three areas of navigation : R=1, R=2 and 
R=3. The verified ship has been distinguished by lack of a lar- 
ge superstructure on the main deck. Only a small deckhouse 
accomodating ship's wheelhouse, has been provided on its bow. 
Calculation results obtained by means of the algorithm in que- 
stion deals solely with hull structure mass, without outfitting, 
therefore the verification serves only as an estimation of pre- 
diction credibility. 

Tab. 1. Results of verifying calculations 


Verified 
quantities 
Mass 
centre 


Values estimated by using the algorithm 


Navigation area | Navigation area | Navigation area 
R=1 R=2 R=3 


(MC) MC MC MC 
height a height a 3 height a height 
[m] [m] [m] [m] 
168.6 


The verifying calculations were performed 
by using Excel calculation sheet. 


RECAPITULATION 


> The presented algorithm for estimating hull structure mass 
of inland navigation ships can be useful in the preliminary 
design stage of inland navigation ships intended for carry- 
ing general cargo, dry and liquid bulk cargoes. In the case 
of ships of different construction, e.g. those for carrying 
heavy roll-on — roll-off cargo, fitted with heavy decks, re- 
sults obtained by means of the algorithm may significantly 
differ from real masses as structural designing such ships is 
carried out in a different way. 


© The performed verifying calculations showed convergent 
results of estimated masses, and, under the assumption that 
the SINE 207 ship construction satisfies the requirements 
for the navigation area R = 2 , the obtained estimation can 
be deemed very correct. 


> The large difference of estimation of mass centre height 
results from not accounted for ship's outfitting. 
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Miscellanea 


Founded in 1994 by the State Committee for Scientific 
Research (KBN), Academic Computer Centre (CI TASK) 
in Gdansk is an inter-university unit that manages one of the 
biggest and most modern metropolitan area networks (MAN) 
in Poland. 

The TASK network covers the territory of the whole so- 
called Tri-City, i.e. Gdansk-Sopot-Gdynia. It connects 70 
LAN networks of various research institutes in which over 
6000 computers of PC class, workstations, and servers are 
installed. The network has about 16000 users (excluding stu- 
dents). The TASK network is connected to the national net- 
work thus enabling remote access to the Centre resources, 
the usage of all networks services, multimedial transmissions, 
interactive work, teaching and learning through the network. 
The Centre co-operates closely with four other supercompu- 
ter centres in Poland, and participates in creation of national 
metacomputing system. It also has a link to the TEN-155 
panEuropean research network which provides extension of 
the service to the USA and other overseas regions. 


A concise characteristics of the Centre 
main departments is as follows : 


The Network Group serves as the design and supervi- 
sing centre of the whole Tri-City network. Now the network 
operates in ATM and FDDI technologies. It consists of three 
FDDI rings, each with five to eight FDDI nodes connected 
with links operating in ATM technology. The group main- 
tains all the network services and servers, monitors secure 
network operation, and gathers network statistics, both 
within TASK, cooperating MANSs, and worldwide transmis- 
sions. Network supervising is highly automated with such 
tools as SUN Net Manager, Cisco Works, and Cisco Net- 
-Flow Analyzer. Lately, the Network Group has offered cre- 
ation of the Virtual LAN facility. A pilot VLAN configura- 
tion has been established for the University of Gdansk. 


The Academic Computer Centre 
in Gdansk (CI TASK) 


The Supercomputers Group maintains and manages 
all the computers (IBM, SGI and SUN) installed in CI TASK 
equipped with a high-capacity archivisation system (ATL, 
HP, and EXABYTE). 34 profesional software packages 
make it possible to perform large scale calculations in va- 
rious fields of science and technology. The group also deli- 
vers expert advice in software application and security of 
servers. Its latest attainment is the SGI cluster project that 
resulted in effective supercomputer workload. 


The Inter-Disciplinary Group of Mathematical Mo- 
delling was established as a group for mathematical mo- 
delling of theoretical and applied research. Its main fields 
of interest include: Quantum Chemistry, Plasma Physics, 
Atomic Physics, Biomolecular Modelling, Mechanics, 
Medicine, and Electrodynamics. The Group publishes pa- 
pers of its interests in the TASK Quarterly Scientific Bulle- 
tin. The journal makes it possible to present papers and 
exchange of views on applied numerical methods for sol- 
ving a variety of problems in science and engineering by 
using high performance computers. Contributors are invi- 
ted to send their proposals to : 


quarterly@task.gda.pl 


In 1999, in recognition of the Centre’s achievements in 
metacomputing technologies based on SGI servers, a spe- 
cial prize was awarded to it by the KBN (The State Com- 
mittee for Scientific Research). 


Notation : 
ATM - Asynchronous Transfer Mode 
FDDI - Fiber Distributed Data Interface 
LAN - Local Area Network 
MAN - Metropolitan Area Network 
VLAN - Virtual Local Area Network 


ATL, EXABYTE, HP, IBM, SGI, SUN - firm names 
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ECODOCK 


The Faculty of Ocean Engineering and Ship Technolo- 
gy, Gdansk University of Technology, having completed the 
Baltecologicalship project, realizes another project within 
the frame of EUREKA European projects. It is E!2968 pro- 
ject whose full name is : 


Environment Friendly Floating Dock 


The program's aim is to perform research and analyses 
to determine probable directions and intensity of future ne- 
eds for environment friendly docks. The measurements and 
systematic investigations of pollution emitted during ship 
repair process will be carried out in ship repair yards. The 
obtained data and assessment of needs will allow to prepare 
technical requirements, initial design and assumptions for 
the construction process of environment -friendly and cost- 
-efficient ecological floating dock of a new type. 

The project will include the in-depth analytical work 
dedicated to the reduction of the level of hazard for the sea, 
coastal area, ports and adjacent estuaries of rivers, assuming 
application of modern technical solutions of proecological 
components and systems, in order to decrease the level of 
harmful emission and amount of waste materials. The expec- 
ted result is to establish design requirements for technically- 
-advanced, environment-friendly and cost-effective floating 
dock of a new type able to safely serve for all kinds of do- 
cking operations and ship repair work, particularly for the 
ECOLOGICAL VESSEL. 

The realization of the project has been divided into seve- 
ral tasks to be realized by co-operating scientific centres and 
enterprises, namely : 


Gdansk University of Technology : 


> Assessment of the present and future market demand for 
ship repair work in Europe versus the situation on the 
world market, determination of future demand for ship 
repair potential ( overhauls and conversions) in the area 
of Baltic Sea. 

> Analysis of quantitative state, condition and potential of 
ship repair industry in the Baltic Sea area. 

> Design assumptions for the ecolgical floating dock and 

an adapted dock, and their developing study. 

Study and compendium of legal, organizational and tech- 

nical regulations concerning floating docks, docking 

operations, as well as building the docks and their con- 

versions with taking into account environment protec- 

tion aspects. 

Development concepts of a new generation of classifi- 

cation rules for construction and operation of floating 

docks. 

Characteristics and classification of typical procedures 

and processes of repair , modernization and conversion 

of ships. 

Proecological operational concepts of service procedu- 

res for a floating dock and realized processes of ship re- 

pair and conversion. 

Analysis of technical solutions and development of a new 

generation of floating docks and their systems with ta- 

king into account operational reliability and safety as well 

as ecological cleanness (design, building, operation, over- 

hauls). 

Analysis of ship docking-in-and-out processes with ta- 

king into account selected critical damage states of ships 


and ecological hazards resulting from them (hydrome- 
chanics, strength, engineering processes). 

> Current and final verification of designs of a new ecolo- 
gical dock and a dock adapted to the CLEAN class. 


Warsaw University of Technology : 


Studies on and compendium of technical and physical 
environment protection problems associated with building 
and operation of floating docks : 
© ecological criteria for floating docks 
= identification and qualification of main ecological pro- 
blems associated with building and operating a classical 
floating dock and CLEAN-class dock 
identification of hazards to the environment from the side 
of ship repair processes 
identification of qualitative and quantitative structure of 
production remains and waste materials from ship repair 
yards, and principles for their management 
procedures for minimizing noxious emissions and waste 
materials in service conditions of ship repair yard 
problems and methods of monitoring the engineering and 
operational processes during ship overhaul and conser- 
vation work as well as during building floating docks, in 
the environment protection aspect 
proecological concepts of procedures and methods for 
floating dock operation, realized during ship repair and 
conservation work 
recommendations concerning devices, systems and engi- 
neering processes for proecological floating docks. 


SINUS Design Office Co : 


> Analysis of technical and ecological shortcomings of a se- 
lected existing dock, choice of the dock as well as a scope 
and structure of proecological modernization work. 

> Elaboration of design assumptions for proecological mo- 
dernization of the selected floating dock. 

> Taking part in elaborating design assumptions for a new 
generation of proecological floating docks of CLEAN 
class. 

> Elaboration of preliminary design of a proecological flo- 
ating dock. 

> Elaboration of preliminary project of proecological mo- 
dernization of the selected existing floating dock. 

> Prediction of building cost at selected technological con- 
ditions. 


Maritim Shipyard, Gdansk : 


* Taking part in elaborating proecological design assump- 
tions of a floating dock intended for own needs. 

* Taking part in monitoring the environment (water, land, 
air) in the area of its own hull building centre. 

* Taking part in monitoring the environment in the area of 
docking centre of NAUTA Shipyard or Gdynia Naval 
Shipyard. 


Schiffsmakler GmbH, Rostock : 
#* Assessment of German ship repair market. 
INFERT GmbH, Rostock : 


«* A study of industrial environment protection systems 
in the range of ship docking processes. 
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ABSTRACT 


In the paper were presented examples of use of some intelligence tools such as a neural network, expert 

system and relational database to ship design. The neural network of back-propagation of errors was 

applied to select required power of ship main propulsion system on the basis of ship main parameters. 

Results obtained by using the network were compared with resulting values for similar ships found in 

Access database application. To aid design of the main propulsion system and ship power plant automation 

fuzzy logic was applied as an element of Case Based Reasoning (CBR) method in Exsys expert system as 
well as a few methods for selection of similar ships, elaborated by the authors. 


Keywords : artificial intelligence, expert systems, neural networks, relational database, 
Case-Based Reasoning method, aided ship design 


INTRODUCTION 


In the subject-matter literature examples of application of 
artificial intelligence to ship design can be found, especially to 
its preliminary phase when ship’s main parameters are selec- 
ted on the basis of shipowner’s design assumptions. In oder to 
ensure optimum main dimensions for a ship during its desig- 
ning the approach which consists in finding ships of similar 
characteristics and modyfing the selected design solutions, is 
often applied. 

To use information dealing with earlier elaborated similar 
designs is possible both by means of expert systems on the 
basis of the Case Based Reasoning method which makes desig- 
ning a ship of high effectiveness faster and easier, or by means 
of neural networks which can be teached on the basis of repre- 
sentative examples and results achieved from other sources (e.g. 
from ship service). Thus the information processing characte- 
ristic for traditional expert systems can be deemed comple- 
mentary to the dispersed parallel processing typical for neural 
networks. 

In this work both the classical artificial intelligence tools 
together with a relational database were used to exemplify the- 
ir application to aiding ship design in the following range : 


* Selection of power output of ship’s main engine (ME) on 
the basis of overall ship’s parameters (mainly its dimen- 
sions) by using similarity calculation methods applied in 
database, based on the Case Based Reasoning (CBR) ap- 
proach, and then possible verification of results with the 
use of a neural network. 

* Aiding design of ship main propulsion system (MPS) by 
means of selection of similar ships with the use of fuzzy 
logic method on the basis of such main propulsion system 
parameters as ME power output and speed, as well as simi- 
larity of the remaining parameters achieved from database 
application. 

* Aiding selection of some parameters of ship power plant 
automation system on the basis of a domain model or exis- 
ting designs of similar ships, with the use of fuzzy logic 
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method embedded in an expert system, in cooperation with 
database application. 


In order to select ME power output on the basis of general 
ship parameters by using similarity calculation methods in da- 
tabase application it was made it possible to verify obtained 
results by means of a neural network which uses the error back- 
-propagation method. And, for similarity calculations and se- 
lection of similar ships within the scope of design of ship main 
propulsion and power plant automation systems the fuzzy lo- 
gic method [9] embedded in Exsys system was used, apart from 
the methods provided in database application [6]. 


SELECTION OF MAIN ENGINE POWER 
OUTPUT BY USING CBR METHOD 
AND NEURAL NETWORK 


To select main engine power output a neural network was 
used apart from the ship similarity calculation methods provi- 
ded in database application. The below presented research re- 
sults are based on the set of 222 ships built by Polish ship- 
yards. The set is composed of very different ships. Each of 
them is characterized by the following parameters : 


Dwt — deadweight, L — overall length, B — breadth, 
D — draught, V — speed, as well as ME rated output. 


In the performed investigations a relationship between ME 
power output and the remaining parameters was searched for. 
A fragment of the ship database containing values of the para- 
meters in question is examplified in Tab.1. 

The ME power output values searched out of the ship data- 
base, most similar to the power chosen according to particular 
similarity calculation methods, were compared with the values 
obtained from the neural network. 

In neural network applications the one-directional networks 
and the teaching method with teacher’s assistance are used most 
often. They were also applied to the investigations on aiding 
ship design. A classical method of teaching the multi-layer, 
one-directional network is the back-propagation algorithm, the 


Tab. 1. An example fragment of the ship database 


Ship] Dwt L v 
No. It] [m] [kn] 


ME rated 
output 
kW 


most often used in technical applications aimed at modelling 
unknown processes. It is the most widely known and used al- 
gorithm for the neural networks of nonlinear output function. 
The algorithm’s essence consists in the reverse direction of 
correcting the weights (teaching the network) : beginning from 
the initial layer to the first hidden layer preceding it, and fur- 
ther up to the first layer. Error measure is a function of net- 
work’s weights. Teaching the network consists in an adaptive 
correction ofall weights in such a way as to obtain a minimum 
of the measure. 

The error back-propagation method has contributed to broa- 
dening the application range of neural networks a.o. to ship 
design. 

In the calculations in question a two-layer network of a con- 
tinuous unipolar activation function and the classical algorithm 
of back-propagation of errors of weight changes was used. The 
set of ships was split in two subsets : teaching and testing. To 
the testing set 25% of ships were chosen in random. A fragment 
of the teaching set is shown in Tab.1. The ships of numbers not 
shown in col.1 of the table were assigned to the testing set. 


All values of the ship parameters were normalized in ad- 
vance to obtain the values within the interval [0,1]. In this case 
one calculation cycle comprised putting-in parameters of all 
ships of the teaching set to the network, successively. The net- 
work teaching was completed when the cycle rms error Em 
achieved a value smaller than the assumed one. The error con- 
cerned the difference between the real ship ME power output 
and that calculated by the network for the same ship. 

Convergence of the teaching process with the use of the 
teaching set of ships, expressed by the relationship between 
Em and the number of calculation cycles Lec, is presented in 
Fig.1. The results were obtained for the network of the follo- 
wing features : six inputs, 25 neurons in the hidden layer, one 
output neuron, the teaching factor n = 0.5, the activation func- 
tion parameter B = 1.5. 


0.035 
0.030 
0.025 f 

E 0.020 
0.015 
0.010 
0.005 


100 
Lee 


150 200 250 


Fig. 1. Teaching process convergence of the applied neural network 


On completion of teaching the network by means of the 
teaching set of ships the calulations based on the obtained net- 
work’s weights were performed with the use of parameters of 
ships from the testing set. The calculated ME power outputs 
are compared with the real ones in Fig.2. In the figure are also 
presented the results obtained by means of the multi-dimen- 
sional regression method with the use of 3° order polynomial 
model [4]. 


In the testing calculations the relative rms error e 
expressed by the formula (1) was also calculated : 


25 000 

20 000 
= Oreal ME power output values [kw] 
= E ME power output values [kW] from regression method 
fal E ME power output values [kW] from the neural network 
n 
3 15000 4 
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3 
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a 
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Ship numbers 


Fig. 2. Comparison of the results obtained for the testing set of ships 
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where : 


M; , Mo — real ME power output 
and that calculated by the network, respectively. 


The following values of the error were obtained : 


e1 = 25.61% for the regression method 
and e2 = 23.13% for the network. 


On the basis of the results of the presented research 
on the neural network use the following conclusions 
can be offerred : 


*# the results obtained from the neural network and those from 
the regression method are similar, respectively 

* the large discrepancy of the results for the testing set (large 
values of the error e) results from the large variety of the 
considered ships 

* for the gradient teaching algorithm of the constant teaching 
factor n a fast convergence in the initial phase and a very 
slow one in the further phase of calculations, is characteri- 
stic (Fig. 1). 


For practical applications more effective teaching algori- 
thms and methods of changing the factor n in the course of 
calculations, should be used. 

In the literature sources various modifications of the error 
back-propagation algorithm, aimed at acceleration of the algo- 
rithm’s convergence, have been proposed [10]. 

On the basis of the following ship design parameters : dis- 
placement, overall length, breadth, draught, and speed the si- 
milar ships were selected and their ME power output values 
were compared with those designed as well as with relevant 
values obtained from the neural network [4]. Calculation re- 
sults dealing with the ships built by Polish shipyards, both ob- 
tained by applying the similarity calculation methods embed- 
ded in the database and those achieved from the neural net- 
work are presented in Tab. 2. 


Tab. 2. ME power output values of similar ships, obtained by using different 
similarity calculation methods, and by applying the neural network 


ME ME power output of similar ship 


power calculated by using the following similarity 
output of methods 


designed based on 
ship function 


[kW] with lower 
limit 


from the 
neural 
network 


based on | based on | based on 
Gaussian |trapezoidal| triangular 
function | function | function 


Cino | oso | soso [12960 [isos [9 


In many cases the designer can assume values of certain 
parameters with a lower or higher tolerance; it specially con- 
cerns the upper limit for a required parameter value as another 
criteria may be decisive, e.g. preference for a given supplier, or 
more favourable terms of delivery. For other parameters to lo- 
wer a limiting value of a designed parameter may be recom- 
mended e.g. with a view of price criterion. From this point of 
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view choice of a form of similarity function is important. In 
the subject-matter literature mainly symmetrical similarity func- 
tions are presented, e.g. the symmetrical similarity based on 
the theory of sets, or the above mentioned symmetrical simila- 
rity with a lower limit [9]. 

These authors have enlarged the set of similarity functions 
by the following ones : triangular, trapezoidal and Gaussian 
[6]. On the one hand it offers, for designer, a greater flexibility 
in determining similarity, on the other hand it forces him to 
select a function complying with requirements and type of an 
analyzed design parameter of automation system. 

The similarity function with lower limit makes it possi- 
ble to limit the range of calculations of similarity ofan investi- 
gated parameter by establishing its lower limit. In this case the 
similarity value linearly increases from zero, at the lower li- 
mit, up to one. 

The trapezoidal and triangular similarity functions make 
it possible to limit the range of similarity calculations by esta- 
blishing two limits : lower and upper, and in the case of trape- 
zoidal function — also the deviations : lower and upper, allo- 
wing to apply some tolerance margin for a given parameter. 

The Gaussian similarity function treats the similar cases 
with the highest tolerance, not rejecting even the least similar 
ones. 

From the presented example it results that the ME power 
output values for the most similar ships, obtained by using par- 
ticular similarity calculation methods, are not always close to 
those of the designed ships. It results from that the similar ships 
were selected on the basis of total similarities ofall input para- 
meters. Therefore it is very important to appropriately esta- 
blish values of weights of the parameters as well as limiting 
values for investigated ranges and their deviations. 

Summing up the performed investigations one can state that 
the best results are obtained by using the calculation method 
based on the Gaussian function, and in some cases — the me- 
thod based on the trapezoidal function. 


The differences of the similarities obtained by means 
of particular methods may result from the following causes : 


Æ a too low number of input parameters of crucial importan- 
ce for choice of ME power output, has been taken for the 
comparisons 

æ% a very diverse structure of the investigated set of ships in 
the database (different types, tasks, gabarites) 

# a too low number of the analyzed ships 

æ% atoo small set of ships in the database, which is especially 
important for teaching the neural network. 


FUZZY LOGIC AIDED DESIGN OF SHIP 
MAIN PROPULSION SYSTEM 


The aiding of design of ship main propulsion system (MPS) 
is carried out by using the probability calculation methods con- 
tained in database application, as well as fuzzy logic method 
embedded in expert system. 

Fuzzy logic can be applied simultaneously to any number 
of parameters. It makes it possible to pass from numerical qu- 
antities to linguistic ones by which logical reasoning can be 
easily performed [8]. The process this way becomes indepen- 
dent of a scale of numerical values of considered parameters. 
A general schematic diagram of the discussed method of pro- 
bability determination is presented in Fig.3. 


Resulting value of numerical measure 
of the similarity Mp depends on : 


+ assumed attribution functions by which input values are 
fuzzified 


+ definition of rules by which reasoning is carried out 
+ mode of sharpening realization. 


R4 Lı 
—|Fuzzification 

R3 Ly 
—|Fuzzification 

R, L 
—".|Fuzzification-— 


Fig. 3. Schematic diagram of the fuzzy set method 


Logical reasoning |+| Sharpening > 


where : 

R,, Ro,..R, — differences of numerical values of conside- 
red input parameters 

Li, L,,.. La — sets of linguistic definitions associated with 
assumed fuzzy sets defined on the basis of 
the values Rj, Ro,.. Ry, respectively 

P, — fuzzified probability 

M, — resulting numerical similarity measure. 


The MPS design aiding is carried out in the following way : 


> in the database application similarity of MPS is calculated 
on the basis of such non-numerical parameters as : type of 
ME and its producer, as well as the following numerical 
parameters : number of MEs and number of propellers. 

> the weighted sum of similarities of the parameters is intro- 
duced, together with the data on ME power output and 
speed, into the expert system Exsys in which the parame- 
ters are fuzzified. 


Assessment of the calculation methods of MPS similarity 
of ships in the expert system was carried out on the basis of the 
data on the selected MPS parameters given in Tab.3, as well as 
the values of the MPS design parameters contained in Tab.4. 


In Tab.4 values of all MPS parameters 
of a designed ship are recorded. 


As far as the MPS is concerned, determination of similarity 
of ships by using fuzzy logic is based on the following nume- 
rical parameters : 

the absolute value of difference of ME power output of 

a designed ship and built one, Rm 
@ the absolute value of difference of ME speed of a designed 

ship and built one, R, 

@ the non-numerical similarity of MPS (calculated in the da- 
tabase) attributed to built ship, P, , taking values from the 
interval [0,1]. 

Calculation of fuzzified similarities 
is realized in two phases : 


> the fuzzified similarity of MPS, P}, is calculated on the ba- 
sis of the fuzzified parameters Rm and Rn 

> the resulting fuzzified similarity of MPS, P,, , is calculated 
on the basis of the numerical similarity of MPS and fuzzi- 
fied similarity P, . 


Ships similar to a designed ship are selected on the basis of 
the maximum resulting fuzzy similarity. To realize the fuzzifi- 
cation process of values of the parameters the attribution func- 
tions shown in Fig.4, were applied. 

The example logical relationships assumed in rules of the 
expert system, dealing with the numerical fuzzified similarity 
of MPS are given in Tab.5, and the resulting fuzzified similari- 
ty of MPS — in Tab.6. 


Tab. 3. Data of MPS parameters selected out of the database 


Main engine (ME) Number 


of 


of units | Propellers 
6ZB40/48 | SULZER 500 
12V40/54A | M.AN. 450 
16V32D [WÄRTSILÄ 800 


IS35MC | B&W 1 
5S600MC | B&W 
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7. | si3smc_ | B&W 
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6RTASS | SULZER 
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Tab. 4. Values of MPS design parameters 
Main engine (ME Number 


Number Output | Speed of 
propellers 
Designed HCP 

sup Lt [TAT | sun [1722 


19. | SBA8M528 
. | SBA8M528 
SBA8M528 
K9Z60/105E 
6K62EF 


DEUT. 
DEUT. 
DEUT. 
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Ship’s 
symbol 


H (Rm) 


[Rm] 


H (Rn) 


[Ra] 


0 25 50 75 1 Po 
Fig. 4. Attribution functions for MPS parameters 
U-value of membership function determining fuzzy set 


The example rules resulting 
from Tab.5 and 6 can be written as follows : 


If (Ry = zero) and (Ry = zero) then (P; = large) 
If (P; = large) and (Po = small) then (Pw = mean). 


In the Exsys system the degrees of attribution to fuzzy sets 
are taken as the so called confidence values. Reasoning by using 
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Tab. 5. Logical relationships in 1“ phase of reasoning 


| Rm | Zeo | Small _ | lage _ | 


| P, [large | large | mean | large | mean | small | mean | small | small | 


Tab. 6. Logical relationships in znd phase of reasoning 


| Po [small | mean | large | small | mean | large | small | mean | large | 
| Ps [mean | large | large | small | mean | large | small | small | mean | 


appropriate rules triggers calculating confidence values con- 
cerning conclusions from the rules, on the basis of confidence 
values for premises. On their basis the sharpening is realized 
in result of which a single numerical value is achieved : 


Z = WC, + W2C2 + W3C3 
where : 


Z — numerical value of MPS similarity (conta- 
ined within the interval [0, 1] 


C4, C2, C3 — confidence values for the values : „large”, 
„mean”, „small”, respectively, achieved 
from the reasoning process 

W1, W2, W3 — weights having non-negative values, 


W,+ W2+ W3=1. 


The application of fuzzy logic was tested on several exam- 
ples (the design data P1+ P5), and their results dealing with 
similarity calculations and selection of similar ships are shown 
in Tab.7. 


Tab. 7. Calulation results of MPS similarity obtained by using Exsys system 


MPS 
speed of 
similar 

ship 


MPS 
Similarity | power of 
similar 

snip p 


Number 
Example _of 
symbol 


Nos. of 
similar 
ships 


[11400] 110 | 20 [0.6286 [ 10800[ 118 | 7 | 
| 6600 | 150 | 1 | 08 | 6650] 154 | 4 | 
| P4 [i000] 120 | 38 | 0.6286] 13050] 124 | 63 | 


In Tab.8 the calculated similarities and selected similar ships 
(examples P1+ P3) selected with the use of fuzzy logic are com- 
pared with those obtained by using the database application. 


| 
Ea AA 
| P | 


Tab. 8. Results of selection of similar ships 
from database and by using Exsys software 


MPS 
power of 
similar 


MPS 
speed of 
similar 


MPS MPS 
power as | speed as 
designed | designed 


PW) | rpm | ikw] | Irpm] | 


— 
6 600 6 650 applied 
methods 


Similar 
ship 
number 


Applied 
method 


Example 
symbol 


93,117, Database 
SS ae 
P3 17 000 500 84 methods 


AIDING SELECTION PROCESS OF SOME 
PARAMETERS OF SHIP POWER PLANT 
AUTOMATION SYSTEM 


Apart from aiding MPS design, fuzzy logic was also ap- 
plied to aiding selection process of some parameters of ship 
power plant automation system.In the case in question the si- 
milarity was related to the features characterizing power plants 
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of built ships as it was assumed that solutions dealing with 
automation system depended on some features of ship power 
plant. Because of the large number of accounted-for features, 
similarity of ships was determined by using some groups of 
the features. Full set of the considered features was split into 
subsets related to : 


> the whole ship (general similarity determined by a type and 
size of ship) - Tab.9 

> the main propulsion system (MPS) (similarity determined 
by a type of propulsion system and its main parameters) 
Tab. 10 

> selected ship systems (similarity determined by system's 
function and its design features) - Tab. 10 

> the electric power plant (similarity determined by a type of 
generating sets and their main parameters) - Tab.11 


Tab. 9. Values of general ship parameters 


Classifi 
Automation 
-cation 
class 
mai 


Number of 
Ship 4 Ship ais 
— Ship type a refrigerated 
Eara 
| e 


a 


B191 ship 


1 504 


roll on- l 
B501 at ship 


B683 49 000 
48000 | o | o | 


AE OO 
E p ship 


Results of similarity calculations within the range of the 
subsets are considered as partial similarities whose sum repre- 
sents the total weighted similarity. 

The similarities calculated in the database application were 
transferred to the Exsys expert system to be fuzzified together 
with the parameters whose similarities were determined by 
using fuzzy logic directly. 

In Exsys system the B191 ship was selected as the most 
similar. The partial similarities of this ship calculated with ta- 
king into account the weights, are given in Tab.12. 

The maximum partial similarities together with symbols of 
relevant ships as well as maximum total similarity of ship, i.e. 
sum of partial similarities were transferred from Exsys system 
to the database. 

Partial similarities of the similar ship (B 191), calculated with 
accounting for appropriate weights, are presented in Tab.13. 
The maximum partial similarities dealing with different ships 
are presented in Tab. 14. 

It should be added that when the entire automation design 
project of the selected ship does not satisfy assumed require- 
ments its particular elements may be taken from other built ships, 
selected on the basis of maximum similarities of the systems. 

The results concerning investigation of similar ships, car- 
ried out by using Exsys system, were obtained on the basis of 
Tab. 14 which contains weighted partial similarities of all ships 
stored in the database. 

Aiding the design of ship main propulsion and automation 
systems, based on CBR methods consists in automatic sear- 
ching out of ship database the ships most similar to a designed 
one [5]. In this scope, the designer selects, out of database fields, 
the parameters on the basis of which at first the partial simila- 
rity of MPS, electric power plant and particular ship's systems, 
then the similarity of the whole ship, together with their rele- 
vant weights, will be calculated. The designer may pass over 


Tab. 10. Values of parameters of MPS and selected ship s systems 


r Number ME ME 


Number of 
transmission 
gears 


Number of | Number of 
Number of : 8 
ropellers Type of propeller(s) | valves in valves in 
P fuel system | bilge system 


pieces 


6L70MC | 16200 | 107 | 1 | 1 | _ off fixed pitch 
6L46 6300 | soo | | of controllable pitch 


6RTA 62-R1 | 11400 | 102 | 2 | 1 of fixed pitch 
oRTA76 | 17220 | 104 | 0 | 1 [| _of fixed pitch 


SRTAG2U | 8670 | 102 | 0 [ 1 | _ offixed pitch 
sséoMC | 10200 | 105 | 2 | ot [| _ of fixed pitch 


Designed ship 


| Ipieces} |- | [kW] | [rpmj | [pieces] | [pieces] | ' 
1 
2 
1 
1 
4 
1 
1 
1 


aoe 
=a 
a 
Sa 
| 4 [| 8ZAL 40S | 23040 | sio | 2 | 2 ofcontrollable pitch | 28 | 33 | 
OOo 
= 


16500 | uo | 1 | 1 | of fixed pitch 


Tab. 11. Values of electric power plant parameters 


r GS1 generating set GS2 generating set shaft generator 
Ship symbol 
number_| type | output | speed | number | type | output | speed | p 


| [pieces] | -| Ikw] | rpm] [pieces] | 


Tab. 12. Results of investigation of similarities of ships 


Similari 


Ship of electric Rae 
symbol Total General ot MES power ot ship's 
system slant systems 


| B222 | 0.1992 | 0.2755 | 0910 | 0.2400 | 0.3160 | 


Tab. 13. Partial similarities of similar ship (B 191) 


Kind of similari Weight of parameter Value of similari 


0.343 
system 
Similarity of electric 0.168 
power plant 
Similarity of ship’s 0.260 
systems 
Tab. 14. Maximum partial similarities dealing with different ships 


Kind of similari Value of similari 


B500 


system 


Similarity of electric B222 
power plant 


Similarity of ship 


systems 


B191 


the phase of calculation of weights, in this case they will be 
assumed equal. 

The weighted similarity both partial and total one of the 
whole ship can be calculated by means of a database applica- 
tion program with the use of any of the proposed method ba- 
sing on one of the applied similarity functions ( that with lower 
limit, trapezoidal, triangular or Gaussian one). In the expert 
system using fuzzy logic, for the searched-out ships most si- 
milar to the designed one, are specified values of its partial 
similarities and of the maximum partial similarities found for 
a given parameter. 
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If ambiguous results are obtained from the above mentio- 
ned methods the designer is able to verify them by using an 
error back-propagating neural network, assuming selected da- 
tabase fields as input and output parameters. 


RECAPITULATION 


Q In recent years intensive research aimed at making use of 
developments of artificial intelligence, a.o. expert systems 
and neural networks, in solving different tasks of ship de- 
sign, have been carried out [2], [3], [4]. 

O These authors attempted to apply an error back-propaga- 
ting neural network for selection of power of ship main 
propulsion system (MPS) on the basis of general ship para- 
meters, a.o. ship main dimensions, in order to find a new 
solution in the case when already found one is not satisfac- 
tory for ship designer. 

O The results obtained from the neural network were compa- 
red with the values resulting for the similar ships found in 
the database. It may be concluded that the neural network 
teached in advance on the basis of solutions of existing ships, 
with the use of general ship design data, can be applied to 
verify results achieved by means of the similarity determi- 
nation methods embedded in the database or expert system. 

O Moreover, a fuzzy logic method was used in Exsys expert 
system to aid design of MPS and ship power plant automa- 
tion, as a tool supplementary for CBR method ; this appro- 
ach is an original and novel solution in the area of ship 
designing. 

O The elaborated system together with the database applica- 
tion can serve as an intelligent designer-friendly tool to aid 
design process of ship automation systems. 
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The Technical University of Hamburg - Harburg and 
Schiffbautechnische Gesellschaft organized the successive, 
9th International Symposium on: 


Practical Design of Ships 
and other Floating Structures 


This widely recognized scientific conference was held 
on 12+17 September 2004 in Lübeck - Travemünde (Ger- 
many). It gave an opportunity for international contacts and 
cooperation of experts to stimulate development of design 
and production technology with a view of effectiveness and 
economy as well as safety improvment of ships and other 
floating objects. 


To realize the idea 142 qualified papers 
were presented during 36 topical sessions on : 


> Design methods > Catamarans and pentamarans 

> Resistance > Floating production systems 

> Design loads > Experimental techniques 

> Ultimate load > Production management 

> Manoeuvring > Slamming and sloshing 

> Methodology > Extreme wave loads 

> Trimarans > Hull girder strength 

> Operation > Springing and torsion 

> Dynamic response > Novel ship concepts 

> Design for safety > Lightweight structures 

> Fatigue > Grounding and collision 

> Seakeeping > Production technology 

> Cavitation > High speed monohulls 
< Stability > Reliability analysis 
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Pacific Rim International Conference on Artificial Intelligence 
(PRICAI 2000), Melbourne, August-September 2000. Lecture 
Notes in Artificial Intelligence. Springer Verlag. 2000 

14. Yager R.R., Filev D.P.: Background for fuzzy modelling and 
control (in Polish). WN-T. Warsaw, 1999 

15.Zakarian V.L., Kaiser M.J.: An embedded hybrid neural network 
and expert system in an computer-aided design system. Expert 
Systems with Applications, Vol.15, 1999 

16.Zbroja S.: Case-Based Reasoning in data bases and knowledge 
bases — selected aspects of formal representation of cases (in 
Polish). 24 National Conference on Computer Methods and 
Systems in Scientific Research and Engineering Design. Cracov, 
October, 1999 
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> Padded drives 
> Propellers 
> Steel sandwich 
> Propulsion 


Design for operation 
Vibration and noise 

Marine engineering 
Computer Integrated Design 
and manufacturing 
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In the broad spectrum of topics the greatest number of 
papers dealt with the following : Fatigue (9 papers), Ultima- 
te load, and Methodology (7 papers each), Design methods, 
Resistance, Operation, Design for safety, Stability, and Pad- 
ded drives (6 papers each), Propellers, and Novel ship con- 
cepts (5 and 4 papers, respectively). 

The presented papers were prepared by experts from 16 
European countries as well as Australia, Brazil, People Re- 
public of China, Egypt, India, Japan, Korea, Tatwan and USA. 


Among them were also four Polish authors 
who presented the following problems : 


* Strength test of steel sandwich panel — by J. Kozak 
(Gdansk University of Technology) 

* Numerical simulation of crash and grounding of inland 
waterway transportation barges — by T. Jastrzebski, 
M. Taczala (Technical University of Szczecin), and 
K. Grabowiecki (CIM-MES Project Co., Warsaw) 

* Efficient freight transport on shallow inland waterways — 
— Results of the INBAT R&D project — by T. Jastrzębski 
(Technical University of Szczecin), as a co-author. 


T. Borzecki from Gdansk University of Technology, 
representing Polish scientific circles, took part in the work 
of the International Standing Committee of the ey 


20 POLISH MARITIME RESEARCH, No 1/2005 


Miscellanea 


A valuable monography 
on electric power quality in isolated systems 


The book recently published reflects technical develop- 
ment in the area of application of semiconductor power ele- 
ments and computer techniques to electric power genera- 
tion, transmission and usage onboard sea-going ships. 


The book titled : 


Assessment of electric power quality in ship 
systems fitted with converter subsystems 


was elaborated by Prof. Janusz Mindykowski of Gdynia 
Maritime University, who adressed it to ship designers, bu- 
ilders, shipyard technical supervision service and ship repa- 
ir teams, as well as to ship crews. Obviously it may serve as 
an important aid for scientific and teaching personnel and 
students of technical university faculties dealing with elec- 
tric engineering. 

Ample content of the book based on the author’s know- 
ledge and experience gained from his multi-year research 
and development acitvity, is divided into 10 chapters conta- 
ining the following topics : 


I. Electric power quality in ship electric power systems. 
Character of disturbances occurring in them. 

Il. Specific character of their exploitation. Electric power 
sources, switchboards and consumers onboard ships. 
Requirements dealing with ship electric networks and 
their service conditions. Role of measurements for dia- 
gnostics and control of electric energy processes. 

Ill. Characteristics of converter systems with accounting 
for electronic power converters, and their influence on 
electric power quality on ships. Disturbance sources 
and mechanisms within such systems. 

IV. Measurement and control instruments applicable to 
ships, requirements for them and assessment of their 
accuracy. 

V. Ship electric power plant as an object of measurements. 
Main relationships describing ship three-phase elec- 
tric power network. Analysis of measurement instru- 
ments applied to assessment of electric power quality 
on ships. 

VI. An attempt to defining : what is quality of electric po- 

wer in ship electric networks, by which factors is it 

determined, and on the basis of which premises sho- 
uld it be assessed ? 

Analysis of influence of electric power quality on ope- 

ration safety and economy of ship systems at ship ser- 

vice conditions. 

VII. Electric power quality indices with the view of ship 
electric engineering features. A synthetic index pro- 
posed as a simple multi-criterial assessment method 
of electric power quality. 

IX. Methods of determination of electric power quality in- 
dices by means of measuring instruments, both analo- 
gous and microprocessor ones. Basic mathematical to- 
ols applied to such instruments. Multi-functional mi- 
croprocessor instruments and analyzers of electric po- 
wer quality, elaborated and manufactured by a team of 
Department of Ship Electric Engineering, Gdynia Mari- 


VU. 


time University, under supervision of the author. Des- 
cription of a special analyzer based on own original 
measurement algorithms and software, as well as on 
positive results of own experimental investigations. 
X. Possible ways to improving quality of electric power 
in the systems in question by means of usage of passi- 
ve, active and hybrid harmonic filters, or damping 
transformers and selected solutions of the systems. 
Considerations on possible control of electric power 
quality in accordance with work safety and exploita- 
tion economy criteria for ship engineering systems. 


In that broad presentation of the matters the author took 
into account the up-to-date state of knowledge in the area of 
measurement technique, informatics, electronics and elec- 
tric engineering. The applied mathematical tools seem to be 
appropriate and the way of presentation of many complex 
problems, accessible. The author also indicated directions 
of further research in the area of ,,intelligent” analyzers and 
computerized control systems. 


The book of Prof. Mindykowski is translated from the 
revised Polish edition: „Ocena jakości energii elektrycznej 
w systemach okrętowych z układami przekształtnikowymi”, 
published in 2001 by the same publisher, i.e. Shipbuilding 
& Shipping Ltd, and shortly presented in Polish Maritime 
Research No. 1/2002. 

In comparison with the original, Polish version, to the 
English edition were introduced only small changes consis- 
ting in removal of editorial errors found in the Polish edition 
as well as some subject-matter changes to Ch. 3, 4 and 8, 
associated with amendments (of 2002) of the rules of Polish 
Register of Shipping and other international classification 
institutions. Some cuts resulting from a specifity of transla- 
tion were also made, however the scope of the book remai- 
ned unchanged. 
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A simulation model of energy distribution 
in ship combustion engine 


Nguyen Hoang 
Gdynia Maritime University 


ABSTRACT 


In the paper a model of energy flow and distribution in ship diesel engine cylinder was 
presented. This is a model of discrete parameters, being a continuation of the author 8 rese- 
arch on simulation of energy processing within ship main propulsion engine [2,3]. The 
model in question makes it possible to calculate energy flow values delivered and transfer- 
red during every working cycle of the engine cylinder. Results of application of the model for 
6ZA408S Sulzer engine installed on B672 ship were also attached. The results were compa- 
red with results of test-stand measurements of the engine, performed at different load levels. 


Keywords : computer simulation, ship diesel engines, dynamic processes, energy balance and distribution 


INTRODUCTORY NOTE 


Knowledge of energy flow and distribution within combu- 
stion engine is important for many reasons. It makes it possi- 
ble to estimate and optimize engine’s efficiency as well as to 
diagnostically model its technical state. In the case of analyses 
concerning ship propulsion system it enables to determine va- 
lues of input energy flows to auxiliary systems and waste heat 
utilization devices. 


PHYSICAL MODEL 


Cylinder’s space together with a working medium contai- 
ned in it forms a thermodynamical system which can be de- 
scribed by means of the equation of the il principle of thermo- 
dynamics. The energy delivered to the cylinder in the form of 
fuel flow is processed, due to combustion process, into flows 
of internal work, thermal energy transferred by combustion 
chamber walls, as well as of engine exhaust gas energy. 

The energy transferred by combustion chamber walls is 
absorbed by media contained in the systems of fresh water co- 
oling cylinders, heads and injectors, as well as of lubricating 
and piston cooling oil. Whereas the exhaust gas energy is trans- 
ferred to drive turbine of supercharging system, to produce ste- 
am in waste heat boiler, as well as to the atmosphere and the 
water cooling head and turbine. Internal work of cylinder’s 
working cycle is transferred through the crankshaft-piston sys- 
tem to the propeller, to drive auxiliary devices of the engine, as 
well as consumed during overcoming friction drag within the 
engine. A small part of the combustion energy is transferred 
through the engine’s elements to the environment. 


The following assumptions have been made : 


> Instantaneous temperature, pressure and gas constant of 
a medium take the same values in every point of the space 
over piston. 

© Heat exchange surface area in the cylinder contains only 
three crucial ones : that of lower plate of cylinder head to- 
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gether with valves, of piston head, as well as changeable 
surface area of cylinder liner which is in contact working 
medium. Average thickness of walls of cylinder liner, lo- 
wer plate of cylinder head, piston head and cylinder block, 
as well as material parameters of the walls, are assumed 
constant. 

> Heat flows only in perpendicular direction to particular 

surfaces. 
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Fig. 1. Schematic diagram of energy balance for the thermodynamical 
system "cylinder". Notation : dHg - input enthalpy increment, dHe, - output 
enhtalpy increment, dQ, - heat emitted during fuel combustion, 
dQ, - heat exchanged with walls, pdV - work done by working medium, 
U - working medium internal energy 


MATHEMATICAL MODEL 


Energy balance in the cylinder 
can be written in the following form : 


dU | dQ; _ dQ, _ dv, dH, ee 
dt dt dt 3 dt dt dt 


(1) 


The derivatives appearing in (1) express rate 
of change of the following flows : 


dU ° 
d =U(1) - of internal energy in the cylinder 
a. = =Q (T) - of heat of combustion of fuel injected 
to the cylinder 
i. = Q. (T) - of heat exchanged with the walls 
T 


of internal work performed by wor- 


dV 
p(t) —— = Ni(t) 
dt king medium 


Hi Ha(t) - of input enthalpy 
T 

oe =H,,(1) - of output enthalpy. 
T 


By taking into account that the internal energy of working 
medium is a function of the temperature T and its chemical 
content F 


du _ ðu dT ou dF 
—+— 2) 
dt oT dt oFdt 
where : 


F -ratio of the combusted fuel mass mg, and the mass, mą , 
of the air contained in the mass of working medium m : 


= Ti , m=Mp+tm (3) 
and 
ou 
Cc, =— 
oT 
the equation (1) can be transformed to the form : 
dQ, dQ, dV , dm, 
= eS ee | eee 4 
dio lja d dt d 3 
dt me, dm, dm du dF | 
a a 
dt dt dF dt 


On its basis the instantaneous temperature of working me- 
dium in the cylinder after the calculation step AT , can be deter- 
mined : 


T(t + At) = T(t) + OT xe (5) 
dt 

The equation (5) together with the state equation of wor- 
king medium make it possible to determine course of tempera- 
ture and pressure in the cylinder in function of time or crank 
angle [3]. To this end is also necessary to determine rates of the 
above mentioned energy flows as well as energy parameters of 
working medium in every calculation step of cycle. The wor- 
king medium flows through inlet and outlet valves in normal 
and reverse direction depending on pressure values in the cy- 
linder and manifolds as well as instantaneous state of opening 
of the valves, were determined. Different cases of the equation 
(1) regarding cycle phases realized in the cylinder, were taken 
into account. 

For instance, in the charging phase only inlet valves thro- 
ugh which fresh air flows in the cylinder, are opened. Depen- 
ding on difference of the pressure p in the cylinder and pg in 
inlet manifold, the following flows through inlet valves take 
place : 


> at pg>p- normal flow 


enters the cylinder - of the specific enthalpy of air, Pi in 
the manifold, resulting in changing the chemical content 
of working medium F and its internal energy u in the cy- 
linder : 


dt m 


(F + 1)F dm 
dt 


dt 


dt ol 


dQ 
dt 


dm 
> atpy<p- reverse flow - during which the mass flow (- s z 


flows out - of the specific enthalpy of working nen h, 

in the cylinder. The chemical content of working medium, 
F, and the internal energy u of the medium in the cylinder 
remain unchanged : 


Dror Seo Tag 
dt dt dt 


The equation (4) takes the following form : 


aE ee 
dt me, \ dt Pi dt dt 


When an engine is in a steady working state, it is possible 
to determine an approximate state of working medium after 
each of the working cycle in a given cycle point. Energy balan- 
ce in the engine can be expressed as follows : 


Qa SENi +Q +Q tU 


dm 


(6) 
where : 
Q q -Tate of the energy delivered to cylinder 
N; - internal work power (indicated power) of cylinder 
Q. - flow rate of heat exchanged with walls 
0. - energy rate of gas exhausted from cylinder 
U - internal energy rate of working medium in cylinder. 


In the steady state the internal energy rate U should be equal 
to zero in an ideal model which is practically unavailable. How- 
ever a model in which the value is relatively small, can be de- 
emed correct. 


The instantaneous flow rate of heat transferred through 
walls, Ò> as well as temperature distribution in walls can be 
determined by means of a model of heat exchange process in 
the cylinder — environment system [3]. Heat exchange — thro- 
ugh walls — between working medium and cooling medium is 
an unsteady process. All the parameters : heat exchange surfa- 
ce area, temperature of walls, as well as the convective heat- 
-transfer coefficient are changeable in function of crank an- 
gle. The heat exchange with cylinder walls was modelled by 
means of the finite element method. The cylinder liner was 
divided into n finite elements in which concentrated temperatu- 
res of cylinder liner wall and of cooling water, are determined. 
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Changeable temperature in the direction of wall thickness x 
are derived from the Fourier equation : 


oT, ƏT, m 
da `a? 
with the 3"! kind boundary conditions on the side 
of working medium and cooling water : 
OT, 
a(T -T,,) = -ad ») 
Ox x=0 
(8) 


ay (a -Tu) =A; 2 = 


where : 

Q, Qw — convective heat-transfer coefficient on the 
side of working medium and on the side of 
cooling water, respectively 

as, Às — temperature equalization coefficient and 


thermal conductance of wall, respectively 
Ts, , Ts. — internal and external cylinder liner wall tem- 
perature, respectively 
5s — wall thickness. 


Due to significant temperature and velocity changes of wor- 
king medium contained in cylinder in different phases of wor- 
king cycle, a value of the coefficient œ is calculated in function 
of crank angle, in accordance with the empirical formula [5] : 


-0.2 0.8 -0.53 
alp)=1304 ° plo) To)" 
0.8 (9) 
T 
| CCa +C Elo) -pa) 
IRA 
where : 
@ [deg] - crank angle 
d [m] - cylinder diameter 
p [bar] - working medium pressure in cylinder 
Cn [m/s] - mean piston speed 
Cı=6.18 - for compression, combustion and expansion 


phase, 2.28 - for charge exchange phase 

C2 = 0.00324 [m/s K] - for integral combustion chamber, 
0.00622 [m/s K] - for divided combustion 
chamber 

Vs [m°] - piston displacement 

Tı, pı ,Vı - working medium parameters at the beginning 
of compression 


Pob [bar] - pressure in cylinder in case of external drive. 


The delivered energy flow is determined as the total energy 
contained in the flow of fuel oil delivered to the cylinder, and 
in the flow of charging air: 


Q, =m,W,+m,h, (10) 

The cylinder internal work flow N; is transferred through 

crankshaft — piston system to drive engine shaft . Its magnitu- 

de depends on working medium pressure in the cylinder and 

on engine rotation speed. The latter is determined by means of 

the ship motion model [2] considered as an energy consumer 
model. 


The flow of energy contained in cylinder exhaust gas, Q., , 
is determined as the enthalpy of cylinder exhaust gas flows : 
e e 


- 11 
Q.x = lee Mex (H) 
The flow of heat exchanged with cylinder walls by wor- 
king medium can be considered as a source of the particular 
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heat flows absorbed by cooling media: that of cylinder liner 
cooling water, (O , of cylinder head cooling water, Q.. of 
piston cooling water, Qn , of lubricating oil, Qa as well as of 


the environment, To) : 
ot 


Q, =Qy.+Q, Ont Qat Oa (12) 
If water temperature distribution within cylinder cooling 


jacket is known [3], the flow of heat absorbed by cylinder co- 
oling water, Que , can be determined from the equation : 


Que = m, Cw [Ba = 


where : 


(13) 


Tw,1 , Tw,n — cooling water temperature at inlet to and outlet 
from cylinder jacket, respectively. 


The flow of heat transferred through cylinder head walls 
in steady state can be determined 
by using the following equation : 


Q, =aA,(T-T,,) (14) 
where : 
Ag,Tog — area and temperature of head’s 


lower plate surface, repectively. 


. The flow of heat absorbed by cylinder head cooling water, 
Q » is determined as the sum of the heat transferred by he- 
ad’s lower plate, Qg» and the heat flow given up by exhaust gas 
to valve and exhaust channel walls. The latter is estimated to 
be 10% of the total heat transferred to cooling medium [1,4]. 


Similarly the flow of heat transferred by piston head walls 
in steady state can be determined by using the equation : 


Qu =AA,(T—T,) (15) 
where : 
a - working medium convective-heat coefficient 
A, T; - area and temperature of piston head surface, 


respectively. 


The flow of heat given up to the environment is determined 
as the heat absorbed by the engine outer surface having the 
temperature equal to that of cylinder block. This is expressed 
by the equation : 


Qai = OLA, (T, 


where : 


-T,,) (16) 


OQ;  - ambient air convective-heat coefficient 
Ay ‚Tp - area and temperature of outer surface of cy- 
linder block, respectively. 


The flow of heat absorbed by lubricating oil, Qa can be 
determined as the rest of thermal energy balance (12): 


Qa =Q, “06 Qu- Q, -Qu (7) 


Fig.2 and 3, show the information on distribution of the 
energy flowing through 6ZA40S engine, given by its manu- 
facturer [6]. In Fig.2. is presented the percentage distribution 
ofthe energy of fuel oil delivered to the engine, into : mechani- 
cal energy, heat transferred to cooling media and environment, 
as well as energy contained in exhaust gases. The balance was 
elaborated for nominal conditions of the engine’s operation in 
the tropical zone and it did not account for the heat recovered 
in waste heat boiler. However in Fig.3 shares of the enumera- 
ted flows are shown in function of the engine’s load, related to 
its rated load. Information on the conditions in which the cha- 
racteristics have been obtained, is lacking. 


Heat from exhaust gas 

Exhaust gas losses 28.8%, required for 
Ss heavy fuel 

oil heating 
approx. 2.5% 


Radiation losses 1.0% . 


Heat equivalent of angine power 44.0% 
Total heat RSRS : 
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= 
Cooler water losses 27.2% NS 


Fig. 2. Sankey s energy flow diagram for 6ZA40S engine 
in nominal conditions acc. [6] 
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Fig. 3. Percentage shares of thermal energy flows through 6ZA40S engine 
in function of its load acc. [6] 


CALCULATION RESULTS 


The model of energy flow and distribution, together with 
that of heat exchange process in the cylinder-environment sys- 
tem, was programmed in Borland Pascal code of Delphi pac- 
kage. Calculations of the model were carried out by applying 
changeable integration step with the use of a PC fitted with 
Pentium II processor [2]. 

In Table the results of calculation of the model of energy 
distribution within cylinder are presented for different engine 
load levels. The following flows and their percentage shares in 
the balance : mechanical energy (cylinder indicated power) 
(item. 3), heat transferred from working medium through com- 
bustion chamber walls (item. 4), heat contained in gas exhau- 
sted from cylinder (item. 5), heat transferred through cylinder 
liner and head walls (item. 7 and 8), heat absorbed by lubrica- 
ting oil (item. 10) , as well as heat given up to the environment 
(item. 11). 

The difference between percentage share of the energy de- 
livered to and that transferred from cylinder in steady state was 
shown in item 13 as the relative error of calculations. In reality 
the value has not necessarily to be equal to zero; it equals inter- 
nal energy increment of working medium. This is due to a non- 
linearity of cyclic work of the engine. 

The calculated share of mechanical energy appeared higher 
by abt. 2% than that presented in the Sankey’s diagram ( Fig. 
2) for 100% engine load. Hence it results that the value given 
in the diagram is the percentage share of the effective power 
measured at the engine shaft. 

From Fig. 2. it results that the value given by the engine’s 
manufacturer as the share in the energy balance of the heat 
absorbed by cooling water, is understood as the heat absorbed 
by all cooling systems ( that of fresh water, lubricating oil, and 
charging air); and the share of the exhaust gas energy is calcu- 
lated for exhaust gas behind the turbine. 

On the basis of the calculation results the energy balance 
and its diagram for the modelled engine was elaborated, as 
shown in Fig.4. 


Calculation results of energy distribution in 6ZA40S engine 
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Fig. 4. Sankey ï energy flow diagram for 6ZA40S engine model 
at 100% load in nominal conditions 


On the basis of comparison of the results presented in Fig.3 
and the results of calculation of the model it can be stated that 
the first has been elaborated for constant rated speed of the 
engine. 


The research on dynamics of the processes 
was presented in the author’s paper [3]. 


CONCLUSIONS 


O In the paper the physical and mathematical model of ener- 
gy flow and distribution within the engine were presented, 
which makes it possible to determine instantaneous values 
of particular energy flows in function of engine load. 


O The model was exemplified by means of 6Z40S engine (in- 
stalled as the main engine on B672 ship) and the research 
results were compared with the test-bed measurements of 
the engine. 
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O The condition of formal correctness of the modelled pro- Indices 
cess, i.e. energy balance of the engine, was satisfied in eve- 


: a - ofair 
ry step of calculations at an acceptably low error. b - of cylinder block 
O The model can be used as a tool for investigation of the d - of inlet 
engine operation process, influence of engine parameters = g X ar gas 
. r : y =0 e 
on the process, as well as for diagnosing the engine techni- fb - of combusted fuel oil 
cal state. E 
i - indicated 
O Investigating the model makes it possible to provide new ot - of environment 
recommendations important for laboratory experiments and s - of wall 
w - of cooling water 


conditions of carrying out empirical tests. 


O The model was used for programming a ship propulsion BIBLIOGRAPHY 


system simulator. 1. Benson S.: The thermodynamics and gas dynamics of internal 

NOMENCLATURE combustion engine, Vol. I. Oxford. Clarendon Press. 1982 
2. Nguyen H.: Simulation model of the energy transformation 

a, - temperature equalization coefficient [m2/s] process in medium-speed diesel engine for ship main propulsion 
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h - specific enthalpy [J/kg] diesel engine cylinder - environment system, Polish Maritime 
H - enthalpy [J] Research, No. 1/2002 
m - mass [kg], 4. Teodorczyk A., Rychter T.: Mathematical modelling of piston 
m - mass flow [kg/s] engine working cycle (in Polish). State Scientific Publishing 
N - power [W] House (PWN). Warszawa, 1990 
p - pressure [Pa] 5. Woschni G., Csallner P.: Zur Vorausberechnung des 
Q - heat, energy [J] Brennverlaufes von Ottomotoren bei gednderten 
Q, - heat exchanged with walls of combustion chamber [J] Betriebsbedingungen. Motortechnische Zeitschrift (MTZ) 43/1982 
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W - calorific value [J/kg] 
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‘onference 


Jubilee Symposium 


On 18+19 November 2004 University of Technology. Its scope was very wide and co- 


25th Symposium on Ship Power Plants vere Mie mlowine aes: 
Design, manufacturing and operation processes of ship 


as hele ean power plants and their equipment with accounting for relia- 
As it can be observed the scientific conference gathering bility, safety and environment protection. 
experts working in the domain of ship powering systems 
has had long tradition. It has been born from the idea of 
yearly meetings to exchange knowledge and experience in 
this area among scientific workers of five Polish Baltic Sea 
coast universities, namely : 


It was reflected in 43 presented papers in preparation of 
which the representatives of Gdynia Maritime University with 
11 papers, and Naval University with 10 papers contributed 
to the largest extent. The remaining universities took part in 
the Symposium: Gdansk University of Technology - with 9 


> Gdansk University of Technology papers, Maritime University of Szczecin - with 8 papers, 
> Technical University of Szczecin and Technical University of Szczecin - with 5 papers. 
=?) Gdynia Maritime University The Symposium program was enriched 
> Maritime University of Szczecin by presentation of products of several firms : 
> Naval University of Gdynia. 
ENAMOR - measurement instruments 
Every 5th year each of them serves as the organizer of UNITEST - educational simulators 


MAN -B&W -ship engines 
Alfa Laval - fuel oil purifying devices 
INTERMASZ - filtration membrane devices. J) 


a sucessive symposium in question. The 25th Symposium 
was organized by the Department of Ship Power Plants, Fa- 
Sila of Ocean Engineering and Ship Technology, Gdansk 


GEGGE 
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Miscellanea 


niSail 


This is an educational program 
of Gdansk University of Technology under the heading : 


Enhancement of the University Education 
Process through Sailing Activities 


The basic idea of the project is to make use of the attrac- 
tive features of sailing for social, organizational and scienti- 
fic activation of the University's students and employees. 


The activation having to start the process of stimulating 
students to be interested in sailing will be realized within 
5 activity groups : 


> The Ist group of activities consists of a cycle of optional 
lectures widening education of students. The three-se- 
mester cycle of lectures will contain such topics as : yacht 
building, physical background of sailing as well as mar- 
keting problems in small undertakings such as sailing 
activities. This scope of themes has to provide not only 
technical knowledge but also to demonstrate a role of 
organizational and managing activities in work of high 
qualified personnel. The lectures in question will be pre- 
pared and given by experienced didactic workers. Di- 
dactic materials will be elaborated and published in 
a multimedial form. 


> The 2nd group of activities covers preparation of pro- 
duction and building process of two sailing yachts (of 
about 5 m in length) intended for navigation both in shel- 
tered waters and inland waterways. The building of the 
yachts has to be an attempt to create an equipment base 
for training, sport competition and turism of students. It 
has also to give an opportunity to do student's professio- 
nal practice as well as to gain skill in collective activity, 
coordination of activities, planning, in order to be aware 
of the role of organization and management in engine- 
er's work. 


> The 3rd group of activities covers sailing instruction 
camps (courses). Their aim is to propagate sport and to 
give an opportunity of gaining yachtsman's degrees. One 
of the courses has to prepare a group of about 20 sailing 
instructors — from among students and the University's 
employees — who would be able to propagate sailing by 
taking care of students intersted in getting the skill of 
sailing. The remaining sailing courses have to give an 
occasion to be acquainted with sailing sport for the gro- 
up of 60 students. A term for participation in the sailing 
camps will be obligatory participation in educational part 
of the project ("No free lunch" policy). 


> The 4th form will be organizational activities associated 
with servicing the "small-business" undertakings. The 
problems will be included into the production process of 
yachts as well as medial servicing process of the project. 
Students and teachers will realize the activities as it will 
give them an opportunity to be acquainted with practical 
aspects of realization of engineering undertakings in 
which not only technology itself but also economy, or- 
ganization and management plays an important role. 


> The 5th activity form will be preparation of an internet 
page and multimedial materials about the project. This 
task will be realized by the University's employees and 
students within the frame of practical classes. The stu- 
dents will have an opportunity to apply the newest infor- 
matic technologies to the realized project. 


An additional activity will be a sociological research on 
perception of sailing among students. It is important to 
know : how much sailing is attractive for students, which 
sailing activities are the most popular, whether sailing is con- 
sidered to be important for their future work and studies. 
The research has to answer the questions : which way 
sailing should be included into educational processes of the 
university, and to which extent it would be possible. 


During realization of the project a set of materials, de- 
signs, analyses and lectures making it possible to adapt the 
project to the conditions of another high school, will be ela- 
borated. 


The project's results will be discussed during 
a seminar finalizing the project's realization. 


Teachers from the universities of the North-coast region, 
sailing instructors, sport organizers and representatives of 
local authorities will take part in the seminar. It will make it 
possible to disseminate the educational ideas on which the 
project is based , to transfer the gained experience to other 
circles, to initiate cooperation of the circles as well as to 
assess potential of technical sports for the widening of edu- 
cational offer of the high schools. 


The Gdansk University of Technology serves as a leader 
of the project, and its didactic and scientific workers as well 
as those of the State High Professional School in Elblag and 


Photo: C. Spigarski 
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OPERATION & ECONOMY 


Some aspects of vibration control 


Part II: An optimal active controller 


Part I : „Active and passive correction” was published in No. 4/2004 of this journal 


Krzysztof Kosowski 
Gdansk University of Technology 


ABSTRACT 


The paper presents a theoretical method for determining the optimal correction to be in- 
troduced in a mechanical system. The active control of harmonic vibrations may be achie- 
ved by applying a control unit which ought to reduce the vibration amplitude of the selec- 
ted elements of the system. The proposed method makes it possible to determine the con- 
troller parameters which provide an optimum value of the chosen quality index. This crite- 
rion includes the reduction of weighted amplitudes of the elements on the one hand, and 
minimizes energy of the control signal on the other. The described method is suitable for 
determination of an optimum controller of turbine rotor vibrations caused by bearing oil whip, bearing oil 
whirl or aerodynamic forces. For the case of rotor self-excited vibrations of aerodynamic type the linear 
model of excitations was compared with the neural network method. 


Key words : active control, mechanical vibrations, optimal controller 


PROBLEM DESCRIPTION All of the vectors f, fz, fw, q, qz, qw in the generalized 
scheme are column vectors composed of the corresponding 

Let us consider a mechanical system with n degrees of fre- vectors for a single frequency: 
edom. Its steady harmonic vibrations at the frequency œj; are 


described by the following equation: fı fz, fw, 
f fz fw 
qi =G; fj (28) f= ?| fz=| | fw= 
where : : f 
fi | - vector of amplitudes of harmonic forces (or moments) fk fz, fw, 
acting upon the inertial elements of system qı qZ; qw, 
qi - vector of displacement amplitude of the inertial elements 
of system q= q2 qz= qZ2 qw= qw: 
G; - system dynamic flexibility matrix. : : : 
When taking into consideration vibrations qk qz, qW, 
with & different values of the frequency œ one can write : ee Ae 
=G-f 29) The matrices G, Q, R, Bw, Cw, Bz, Cz 
q are described by the following formulas : 
The active control of the system vibrations is achieved by G 0 0 o 0 0 
applying a controller. The block diagram of the system and its I a= l 
control unit are presented in Fig.11 (for details see Part I of the I0 GĠG = 0 _|0 Q = 0 
G= Q = 
paper). : ae l ; aee 
o o ü 0 0 Qk 
R 0- 0 Bw 0 0 
R0 R - 0 | R-}0 Bw 0 
Fig 11. Block diagram of active control of mechanical vibrations 0 0 © R 0 0 -= Bw 
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Cw 0 0 Bz, 0 0 
Čw=| 0 Owe O | po Be 0 
0 0 Cw 0 0 Bz, 

Cz, 0 + 0 

Gaji SB! 

0 O =: Oy 


It is assumed that external forces fz act on selected inertial 
elements in a way which can be described by the binary matrix 
Bz. The vector qz represents the vibration amplitudes whose 
value has to be minimized. The selection of these amplitudes 
from the vector q is performed by means of the binary matrix 
Cz. In real mechanical systems it is possible to measure the 
amplitudes qw of vibrations of certain elements only. The vec- 
tor qw can be obtained from vector q by multiplying it by the 
binary matrix Cw. Based on the vector qw, the controller (cha- 
racterized by the parameters matrix R) creates the output si- 
gnal fw, by using the transfer matrix R -®. This signal is intro- 
duced to some possible elements of the system, selected by the 
binary matrix Bw. 


The system presented in Fig.11 
can be described by the following set of equations : 


q=G-f (30) 
f = Bz -fz + Bw -fw (31) 
qz=Cz-q (32) 
qw =Cw-q (33) 
fw =R-Q-qw (34) 

k 
R= > [e xI Re x1,,)| (35) 
i=l 
In the case of a PID controller, 
the matrix R takes the following form : 
R=[K,|K,|Kp] (36) 


where : 


Kp, Kj and Kp are matrices of the proportional (P), 
integrating (I) and differentiating (D) action, respectively. 


The optimization problem is to determine a controller ma- 
trix R which minimizes the amplitudes qz of chosen elements, 
taking into account the role of particular elements. This is ob- 
tained by assigning, to each of them, an appropriate coefficient 
from the weight matrix ®. On the other hand, the energy of the 
control signals fw must be taken into account by applying the 
weight matrix A. Thus, the optimization index E may be writ- 
ten as follows : 


E =tr(qz-®-qz+fw-A- fw) 


where : 


(37) 


tr denotes the trace of a matrix, while the superscript asterisk 
of a matrix symbol stands for a transposed conjugate matrix. 
From the set of equations (30+35) 
the following relation can be derived : 


fw =Z-Cw-G-Bz.-fz (38) 


where the matrix Z is defined by the formula : 


Z=R-Q(1-Cw-G-Bw-R-Q)! = 


3 Ske: xI, )R(ej xl,)|2- (39) 
i=l 


k -I 
1-6n-G-nw Se xI )R-Q(e; <ko] 
i=l 
Furthermore, from the equations (30+35,38) 
it is possible to show that : 
q =G(1+Bw-Z-Cw-G)Bz-fz (40) 
The equations (32) and (39) lead directly to the relation : 


qz=Cz-G(I+Bw-Z-Cw-G)Bz-fz (41) 


By substituting fw and qz in the equation (37) for the right 
sides of the equations (38) and (41) respectively, the quality 
index E may be expressed directly as a function of the parame- 
ters of the mechanical system and the controller, the system 
input fz and the weight matrices ® and A: 


E= tr{fz- Bz [(1+G: Cw- Z“ Bw)G: Cz: ®-Cz:G: 
-([+ Bw: Z-Cw-G)+G'Cw-Z'-A-Z: Cw-G]Bz- fz} 


The optimal value of the index E for a controller matrix R 
must occur in a stationary point of the function E(R), i.e. for 
a matrix R, which fulfills the relation : 


OE 
20 
oR 

The condition for all partial derivatives with respect to sin- 


gle coefficients OR[i, j] of the matrix R may be written in the 
form : 


(42) 


(43) 


SRJ j] OZ[k.1] ARG, j] 


It is possible to conclude from the relations (39) and (44), 
that the condition (43) will be always fulfilled for k=1 (which 
refers to the case in which a single frequency is considered ) as 
long as : 


ôE 5 OE OZEN o ay 


oE 
= =0 45 
The formula of the derivative A can be obtained from the 
equation (42), with the aid of fanais given in [1,3]. The re- 


sult may be written in the following form : 
Bw-G-Cz-®-Cz+ 
2 i 
+(Bw-G-Cz-®-Cz-G- Bw + A)Z-Cw 


2E 
OZ 
-G -Bz -fz -fz -Bz - G- Cw (46) 


From the condition (45) for the extremum of function E(Z), 
the matrix Zopt matching the optimal controller matrix Ropt 
can be determined : 


Zop =- (Bw -G-Cz-®-Cz-G-Bw+A)'- 
-Bw-G-Cz-®-Cz-G-Bz- fz -fż -G*. Cw: 
-(Cw- G- Bz: fz-fz-Bz-G-Cw)! (47) 
The matrix Zopt exists if the two matrices inverted in the 


* * 


above given equation, i.e.: (Bw-G-Cz:-®-Cz-G-Bw+A) 
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and (Cw -G -Bz-fz-fz-Bz-G-Cw) are nonsingular. If the 
conditions are fulfilled, the optimal matrix Ropt exists and may 
be easily calculated from the relations (39) and (47) : 


Rope = (l+ Zope’ Cw G- Bw) Zop 2! (48) 


The final values of the controller matrix R should be deter- 
mined by taking into account some additional limits or crite- 
ria, for example stability requirements. 


EXAMPLES OF APPLICATION 
A ship propulsion unit 


opt opt 


An example of a ship propulsion unit equipped with flexi- 
ble couplings, shown in Fig.12, was described in Part I of this 
paper. It was shown that the controller using the angular velo- 
city of the generator as its correction signal offers the largest 
possibilities of reducing torsional vibration amplitude in the 
main coupling, Fig.13. 

The linear dynamic model of this power plant was compa- 
red with experimental data, followed by an analysis of the vi- 
brations due to operation of the engine with one misfiring cy- 
linder [2]. By using the method presented in this paper, the 
optimal matrix R of the parameters of the additional controller 
was estimated, under the assumption that the values of para- 
meters should be contained within a range of practical applica- 
tions. The optimum parameters of the matrix R were determi- 
ned as those fulfilling the criterion (43) for minimizing the qua- 
lity index E in the range of vibration frequencies 12+26 Hz. 

The comparison of the amplitude of the fundamental har- 
monic frequency of the torque acting on the main flexible cou- 
pling is presented in function of engine speed in Fig.13 for the 
cases with and without correction. By using the method de- 
scribed in Part I and II of the paper, it turned out possible to 
change the structure and parameters of the main engine control 
system and to remarkably decrease torque amplitudes (much 
below the hazardous limit). 


GENERATOR} N 


GEAR 


PROPELLER 


MAIN ENGINE 


EFI 


MAIN COUPLING 


p------------------------------------- 


y 


CORRECTION MAIN ENGINE 
CONTROLLER CONTROLLER 


Fig 12. Ship propulsion unit with the correction loop 
for reduction of main coupling vibrations 


Turbine rotor self-excited vibrations 


Currently conducted work (see Part I of this paper) con- 
centrates on active control of rotor vibrations of a 200 MW 
steam turbine by means of pressurized bearings. The aim is to 
investigate the possibilities of applying the pressurized bearings 
to the large output steam turbine rotor to reduce self-excited 
vibration caused by oil whip or oil whirl, as well as aerodyna- 
mic excitations. 


The proposed method is used for : 


> determining the optimum construction of the applied pres- 
surized bearings (i.e. location of the external oil supply, 
Fig. 14a) 
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Fig 13. Comparison of the main coupling torque (current torque 
related to the maximum torque for the case without correction) 


500 


> finding the optimal controller for active control of the con- 
struction parameters of the bearing supports, Fig.14b 

> determining the optimum structure of the active control 
system and the optimum controller parameters for gover- 
ning the rotor vibrations by means of pressurized bearings, 
Fig.15. 


After research is completed, the results 
will be presented in a separate paper. 


Fig 14. Diagrams for determining optimum location 
of external oil supplying (a) and for determining optimum controller 
for active control of bearing supports (b) 


Fig 15. Diagram for determining the optimum structure 
of the active control system for governing the turbine rotor vibrations 


CONCLUSIONS 


Q A theoretical approach for determining the optimal correc- 
tion loop for improving the dynamic behaviour of a linear 
mechanical system was described in the paper. The active 
control of harmonic vibrations was achieved by applying 
a control unit which reduces the vibration amplitude of se- 
lected elements of the system. This method was applied to 
a ship propulsion unit in order to change the structure and 
parameters of the main engine control system, and by use of 
it a remarkable decrease of torque amplitudes was achieved. 


OQ The presented method makes it possible to determine the 
controller parameters which provide the optimum value of 
the chosen quality index. This criterion includes the reduc- 
tion of weighted amplitudes of the elements on the one hand, 
and minimizes energy of the control signal on the other. 

Q The proposed method is suitable for the determination of 
optimum controller of turbine rotor vibrations due to bea- 
ring oil whip, bearing oil whirl or aerodynamic forces. 
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NOMENCTLATURE 

Bw, Bz, Cw, Cz - binary matrices 

e;  - i-th versor of dimension k 

E  - optimization index 

f - vector of force amplitudes 

fw - controller output vector 

fz - vector of external forces 

G - dynamic flexibility matrix 

I - unitary matrix 

j,k, 1,n,1,s - dimensions of matrices and vectors 
K - matrix of stiffness coefficients 

q  - vector of displacement amplitudes 

qw - vector of measured amplitudes 

qz - vector of amplitudes of controlled elements 
R - controller matrix 


®,A - weight matrices 


œ - frequency 
Q - frequency multiplier matrix 
Indices 


P - proportional controller 
I - integrating controller 
D - differentiating controller 
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a C onference — 


21st Scientific Session 
on Marine Technology 


On 27-28 September 2004 Polish shipbuilders met in 
Gdańsk to held their 21* scientific conference which was 
organized by the Polish Society of Naval Architects and 
Marine Engineers KORAB, Faculty of Ocean Engineering 
and Ship Technology of Gdańsk University of Technolo- 
gy (GUT), (which hosted the conference), and North Ship- 
yard Co. 

The Conference started with the session on History of 
Shipbuilding which contained six papers prepared by scien- 
tific workers of the Faculty of Ocean Engineering and Ship 
Technology, GUT, namely : 


* Main Directions of the Development of Ships over the 
Centuries — by K. Rosochowicz, D. Duda 

% A short Pictorial History of Orthodox and Unortho- 

dox Hydrodynamic Marine Propulsors — by J.A. Szan- 

tyr 

Outline of Marine Piston Engines’ Technology Deve- 

lopment — by J. Girtler 

Steam Turbine as the Main Ship Propulsion. More than 

a Hundred Years at Sea — by K. Kosowski 

Gas Turbine : Advanced Marine Propulsion 

by K. Kosowski 

Marine Turbines : Cogeneration and Combined Pro- 

pulsion Systems — by K. Kosowski 


* k*k HH 


The interesting essence of the papers was enriched 
by presentation of many illustrations and drawings. 


39 papers were read during the technical part of 
the Conference, which was divided into 7 topical groups : 


Shipbuilding Techniques (7 papers) 

Ship Design and Hydromechanics (6 papers) 

Ship Hydromechanics and Structural Mechanics 

(7 papers) 

Safety of Ships and Shipping (4 papers) 

Ship Power Plants (6 papers) 

Ship Power Plants and Ship Equipment (4 papers) 
Ship Equipment and Production Management (5 papers). 


+ + > +++ 


Representatives of Gdansk University of Technology 
and Technical University of Szczecin prepared the gre- 
atest number of papers (13 and 11, respectively). Authors 
and co-authors of the! 5 remaining papers came from Gdy- 
nia Maritime University, Polish Register of Shipping, Na- 
val University of Gdynia, Ship Design and Research Cen- 
tre — Gdansk, Wrocław University of Technology, Foun- 
dation for Safety of Navigation and Environment Protec- 
tion, Gdansk Shiprepair Yard. 


tl : C. Spigarski 
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On 12-15 September 2004 it was 30" time 
when the international scientific conference on: 


Internal Combustion Engines 
was held. 


This time the conference venue was Zakopane, a Polish 
health & recreation resort . It was organized by Institute of 
Aeronautics, Warsaw, Research & Development Enterprise 
„PRO-MO” Cracow, and Polish Academy of Sciences. 


The Conference program contained 65 papers 
including 9 ones presented during the plenary session on : 


> Perspectives and trends in ,, Diesel emissions and con- 
trol” — by Thimoty V. Johnson (Corning Incorporated, 
USA) 

> Experimental study on the urea SCR system for medium 
Duty Track — by Tazuki Watanabe, Jin Kusaka (Gradu- 
ate School of Science and Engineering, Waseda Univer- 
sity, Tokyo, Japan) 

> TARDEC- National Automotive Center (NAC) Research 
and Development Overview — by Alexander Sandel 
(TARDEC The National Automotive Center — US Army, 
Waren, USA) 

> Novel pistons for diesel engines — by Jerzy Nykiel (Fe- 
deral Mogul S.A., Gorzyce, Poland) 

> Space vehicles propulsion — some concept in reality and 
fiction — by Cezary Szczepaniak (Lodz University of 
Technology, Poland) 

> Facilitation of HCCI combustion of biogas at moderate 
compression ratios by application of fuel reforming and 
inlet air heating — by M. L. Wyszynski, T. Megaritis, 
J. Karlovsky, D. Yap, S. Peucheret, R. S. Lehrle (The 
University of Birmingham, UK), H. Xu (Yaguar Cars, 
Coventry, UK), S. Golunski (Johnson Matthey Techno- 
logy Centre, Reading, UK) 

> Recent advances in engine air cleaners design and eva- 
luation — by T. Jaroszezyk, Ch. E. Holm (Fleetguard Inc. 
Stoughton, US) and Byron A. Pardue (Fleetguard Inc. 

K Cookeville, US) 


C onference 
KONES 2004 


above enumerated themes, during which 33 elaboration were 
presented (14, 12, 7, respectively). 


reign scientific centres of USA (7), Lithuania (5), Japan (5), 
Great Britain (4), Slovenia (3), Czech Republic (2), Slova- 
kia (2), India (2), Germany (1), Canada (1) and Switzerland 


S 


> Effect of fuel temperature and ambient pressure on a com- 
mon rail rape-seed oil spray — by M. T. Bialkowski, 
T. Pekdemir, D. P. Towers, R. Reuben (Heriot-Watt Uni- 
versity, Edinburgh, Scotland), M. Brautsch (University 
of Applied Sciences, Amberg-Weiden, Germany) and 
G. Elsbett (Elsbett Technologies GmbH, Thalmassing, 
Germany) 

> Effect of composite material on performance and emis- 
sion characteristics of SI. engine - by P. Mahendra Ku- 
mar (National Engineering College, Tamilnadu, India) 
and C.G. Saravanan (Annamalai University, India). 


The remaining papers were read during 3 panel sessions 
devoted to the following themes : 


> Ecology, Combustion, Thermodynamic processes (19 


papers) 
© Fuelling, Measurement, Control (20 papers) 
= Construction, Technology, Operating (17 papers). 


Moreover, had place a poster session divided into the 


Apart from Polish participants 33 representatives of fo- 


(1) took part in the Conference. 


Photo: C. Spigarski 


Mount Giewont (1909 m), a characteristic peak towering over Fane) 


Miscellanea 


IFFM Publishers in Gdansk is a unit of Szewalski Insti- 
tute of Fluid-Flow Machinery (IFFM), Polish Academy of 
Sciences (PAS). 

The publishing house issues the Institute's periodicals 
as well as publications of other bodies, for instance Ther- 
modynamics and Combustion Committee and Committee 
on Energy Problems, PAS, as well as monographies and 
conference proceedings ordered by universities of Gdansk 
sea-coast region. 


Among the Institute's peridicals 
the following should be enumerated : 


@ Transactions of IFFM - the periodical which has been 


IFFM PUBLISHERS 


issued two times a year since 1960, since 1997 only in 
English 

@ Scientific Installments of IFFM — (Studia i Materiały, in 
Polish) — publishing series which presents selected artic- 
les in the form of separate dissertations 


@ Annual Report (Przegląd Prac) — (both in Polish and En- 
glish) — the yearly report on achievements of IFFM scien- 
tific workers, together with comprehensive bibliography 


% Archives of Thermodynamics, a quarterly in English 

“* Archives of Energetics, a periodical mainly in Polish are 
permanent publications of the above mentioned Com- 
mittees of Polish Academy of Sciences. 
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